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^ Abstract 

V 

^ This article is a summary of four introductory lectures on "Neutrino Experiments," 

_G given at the 2006 TASI summer school. The purposes were to sketch out the present 

questions in neutrino physics and to discuss the experimental challenges in addressing 
them. This article concentrates on specific, illustrative examples rather than providing 
a complete overview of the field of neutrino physics. These lectures were meant to lay 
the ground-work for the talks which followed on specific, selected topics in neutrino 
physics. 
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This article is a summary of four introductory lectures on "Neutrino Experiments," 
given at the 2006 TASI summer school. The purpose was to sketch out the present 
questions in neutrino physics, and discuss the experiments that can address them. The 
ideas were then explored in depth by later lecturers. 

This article begins with an overview of neutrinos in the Standard Model and what 
we know about these particles today. This is followed by a discussion of the direction 
of the field, divided into the three themes identified in the APS Study on the Future of 
Neutrino Physics. 1 This APS study represented the culmination of a year-long effort 
by the neutrino community to come to a consensus on future directions. The report is 
recommended reading for students, along with the accompanying working group white 
papers, especially the Theory Group Whitepaper. 2 

While these lectures used the APS Neutrino Study themes as the core, the emphasis 
here is different from the APS report. The point of a summer school is to teach specific 
ideas rather than provide a perfectly balanced overview of the field. The result is that, 
with apologies, some experiments were necessarily left out of the discussion. Students 
are referred to the Neutrino Oscillation Industry Website 3 for a complete list of all 
neutrino experiments, by category. 

1 Neutrinos As We Knew Them 

Neutrinos are different from the other fermions. Even before the recent evidence of 
neutrino mass, neutrinos were peculiar members of the Standard Model. They are the 
only fermions 

• to carry no electric charge. 

• for which we have no evidence of a right-handed partner. 

• that are defined as massless. 

These ideas are connected by the fact that, unlike other spin 1/2 particles, neutrinos 
can only interact through the weak interaction. 

Even though the Standard Model picture is now demonstrably wrong, this theo- 
retical framework provides a good place to start the discussion. This section begins 
by expanding on the Standard Model picture of the neutrino sketched above. It then 
discusses how neutrinos interact. This is followed by an overview of neutrino sources 
and detectors. 

1.1 Neutrinos in the Standard Model 

Neutrinos are the only Standard Model fermions to interact strictly via the weak in- 
teraction. This proceeds through two types of boson exchange. Exchange of the Z° is 
called the neutral current (NC) interaction. Exchange of the W ± is called the charged 
current (CC) interaction. When a If is emitted, charge conservation at the vertex 
requires that a charged lepton exits the interaction. We know the family of an incom- 
ing neutrino by the charged partner which exits the CC interaction. For example, a 
scattered electron tags a v e interaction, a fi tags a interaction, and a r tags a v T 
interaction. The neutrino always emits the W + and the antineutrino always emits the 
W~ in the CC interaction. In order to conserve charge at the lower vertex, the CC 
interaction is flavor-changing for target quarks. For example, in a neutrino interaction, 
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e en p n up hadrons 

Elastic QuasiElastic Resonance DIS 



Figure 1: Examples of the four types of neutrino interactions which appear throughout this 
discussion and are defined in sec. 1.2. The first two diagrams show an NC and CC interaction, 
respectively. 



if a neutron, n, absorbs a W + , a proton, p, will exit the interaction. The W has 
converted a d quark to a u quark. The first two diagrams shown on Fig. 1 illustrate a 
NC and a CC interaction, respectively. 

In 1989, measurements of the Z° width at LEP 4 and SLD 5 determined that there 
are only three families of light-mass weakly-interacting neutrinos, although we will 
explore this question in more depth in section 3 of these lectures. These are the u e , 
the and the v T . The interactions of the v e and Vn have been shown to be consistent 
with the Standard Model weak interaction. Until recently, there has only been indirect 
evidence for the v T through the decay of the r meson. In July 2000, however, the 
DoNuT Experiment (E872) at Fermilab presented direct evidence for v T interactions. 6 

Within the Standard Model, neutrinos are massless. This assumption is consistent 
with direct experimental observation. It is also an outcome of the feature of "handed- 
ness" associated with neutrinos. To understand handedness, it is simplest to begin by 
discussing "helicity," since for massless particles helicity and handedness are identical. 

For a spin 1/2 Dirac particle, helicity is the projection of a particle's spin (X) along 
its direction of motion p, with operator I] • p. Helicity has two possible states: spin 
aligned opposite the direction of motion (negative, or "left helicity" ) and spin aligned 
along the direction of motion (positive or "right helicity"). If a particle is massive, 
then the sign of the helicity of the particle will be frame dependent. When one boosts 
to a frame where one is moving faster than the particle, the sign of the momentum will 
change but the spin will not, and therefore the helicity will flip. For massless particles, 
which must travel at the speed of light, one cannot boost to a frame where helicity 
changes sign. 

Handedness (or chirality) is the Lorentz invariant (i.e., frame-independent) ana- 
logue of helicity for both massless and massive Dirac particles. There are two states: 
"left handed" (LH) and "right handed" (RH). For the case of massless particles, in- 
cluding Standard Model neutrinos, helicity and handedness are identical. A massless 
fermion is either purely LH or RH, and, in principle, can appear in either state. Mas- 
sive particles have both RH and LH components. A helicity eigenstate for a massive 
particle is a combination of handedness states. It is only in the high energy limit, 
where particles are effectively massless, that handedness and helicity coincide for mas- 
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sive fermions. Nevertheless, people tend to use the terms "helicity" and "handedness" 
interchangeably. Unlike the electromagnetic and strong interactions, the weak interac- 
tion involving neutrinos has a definite preferred handedness. 

In 1956, it was shown that neutrinos are LH and outgoing antineutrinos are RH. 7 
This effect is called "parity violation." If neutrinos respected parity, then an equal 
number of LH and RH neutrinos should have been produced in the 1956 experiment. 
The fact that all neutrinos are LH and all antineutrinos are RH means that, unlike all 
of the other fermions in the Standard Model, parity appears to be maximally violated 
for this particle. This is clearly very strange. 

We need a method to enforce parity violation within the weak interaction the- 
ory. To this end, consider a fermion wavefunction, ip, broken up into its LH and RH 
components: 

1p = IpL + TpR- (1) 

We can introduce a projection operator which selects out each component: 



T^L,R = T^L,R- 



(2) 



To force the correct handedness in calculations involving the weak interaction, we can 
require a factor of (l- 7 5 )/2 at every weak vertex involving a neutrino. As a result of 
this factor, which corresponds to the LH projection operator, we often say the charged 
weak interaction (W exchange) is "left handed." 

Note that by approaching the problem this way, RH neutrinos (and LH antineu- 
trinos) could in principle exist but be undetected because they do not interact. They 
will not interact via the electromagnetic interactions because they are neutral, or via 
the strong interaction because they are leptons. RH Dirac neutrinos do not couple to 
the Standard Model W, because this interaction is "left handed," as discussed above. 
Because they are non-interacting, they are called "sterile neutrinos." By definition, 
the Standard Model has no RH neutrino. 

With no RH partner, the neutrino can have no Dirac mass term in the Lagrangian. 
To see this, note that the free-particle Lagrangian for a massive, spin 1/2 particle is 



C = itp^^d^ip — mtpip, 

However, iptp can be rewritten using 

1> L , R = 1/2(1 T 7 5 R 
$ l ,r = 1/2^(1 ± 7 5 ), 



(3) 



(4) 
(5) 



giving 



1 + 7 5 1 - 7 5 



1 + 7 5 | 1 
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1p = 4>LlpR + 1pRi>L- 



(6) 



In other words, an rm/jip ("mass") term in a Lagrangian mixes RH and LH states of 
the fermion. If the fermions have only one handedness (like vs), then the Dirac mass 
term will automatically vanish. In the Standard Model, there is no Dirac mass term 
for neutrinos. 
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Figure 2: Current status of CC cross section measurements in the 1 to 100 GeV range. 
This plot shows cr/E, thus removing the linear energy dependence at high energies. Note 
the low energy cut-off due to the muon mass suppression. Components of the total cross 
section are indicated by the curves. 8 

1.2 Neutrino Interactions 

Neutrino interactions in the Standard Model come in four basic types. Fig. 1 shows 
examples of the four interactions. In Elastic scattering, "what goes is what comes out," 
just like two billiard balls colliding. An example is a NC interaction where the target is 
does not go into an excited state or break up, e.g., v e +n — > v e + n. A more complicated 
example is electron- neutrino scattering from electrons, where the W exchange yields a 
final state which is indistinguishable from the Z exchange on an event-by-event basis, 
so this is categorized as an elastic scatter. Quasi-elastic scattering is, generally, the 
CC analogue to elastic scattering. Exchange of the W causes the incoming lepton and 
the target to change flavors, but the target does not go into an excited state or break 
apart. An example is + n — > (jl + p. Single pion production may be caused by 
either NC or CC interactions. In resonant single pion production, the target becomes 
a A which decays to emit a pion. In coherent scattering, there is little momentum 
exchange with the nucleon and a single pion is produced diffractively in the forward 
direction. The case of NC single 7r° production is particularly important, because 
this forms a background in many neutrino oscillation searches. Finally, DIS, or Deep 
Inelastic Scattering, is the case where there is large 4-momentum exchange, breaking 
the nucleon apart. One can have NC or CC deep inelastic scattering. 

Fig. 2 summarizes the low energy behavior of cr/E for CC events (solid line), as pre- 
dicted by the NUANCE neutrino event generator. 9 The quasi-elastic, single pion and 
deep inelastic contributions are indicated by the broken curves. The data indicate the 
state of the art for neutrino cross section measurements. One can see that if precision 
neutrino studies are to be pursued in the MeV to few GeV range, that more accu- 
rate measurements are essential. The MiniBooNE, 10 SciBooNE, 11 and MINERvA 12 
experiments are expected to improve the situation in the near future. 

Above a few GeV, the total neutrino cross section rises linearly with energy. The 
total cross section is the sum of many partial cross sections: quasi-elastic + single pion 
+ two pions + three pions + etc. As the energy increases, each of these cross sections 
sequentially "turns on" and then becomes constant with energy. Thus the sum, which 
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Figure 3: The ratio of the z/ e to CC cross sections as a function of neutrino energy, showing 
the suppression due to the lepton mass. 

is the total cross section, increases continuously and linearly with E. 

Nevertheless, even at high energies, this interaction is called "weak" for good reason. 
The total cross section for most neutrino scattering experiments is small. For 100 GeV 

interactions with electrons, the cross section is ~ 10~ 40 cm 2 . For 100 GeV 
interactions with nucleons, the cross section is ~ 10~ 36 cm 2 . This is many orders of 
magnitude less than the strong interaction. For example, for pp scattering, the cross 
section is ~ 10 -25 cm 2 . The result is that a 100 GeV neutrino will have a mean free 
path in iron of 3 x 10 9 meters. Thus most neutrinos which hit the Earth travel through 
without interacting. It is only at ultra-high energies that the Earth becomes opaque 
to neutrinos, as discussed in sec. 3.3.2. 

In principle, the interactions of the v e , and v T should be identical ("universal"). 
In practice, the mass differences of the outgoing leptons lead to considerable differences 
in the behavior of the cross sections. In the CC interaction, you must have enough CM 
energy to actually produce the outgoing charged lepton. Just above mass threshold, 
there is very little phase space for producing the lepton, and so production will be 
highly suppressed. The cross section increases in a non-linear manner until well above 
threshold. Consider, for example, a comparison of the v e and CC quasielastic cross 
section on carbon, shown in Fig, 3. At very low energy the CC cross section is 
zero, while the v e cross section is non-zero, because the 105 MeV muon cannot be 
produced. The ratio approaches one at about 1 GeV. A similar effect occurs for the v T 
CC interaction cross sections. The mass of the r is 1.8 GeV, resulting in a cross section 
which is zero below 3.5 GeV and suppressed relative to the total CC scattering cross 
section for v T beam energies beyond 100 GeV. At 100 GeV, which corresponds to a 
center-of-mass energy of \JlME « 14 GeV, there is still a 25% reduction in the total 
CC f T interaction rate compared to v„ due to leptonic mass suppression. 
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Common Terminology 


Reaction 


u pp neutrinos" 


p + p - 


-> 2 H + e~ + v e 


u pep neutrinos" 


p + e~ 


+ p -> 2 H + u e 


" 7 Be neutrinos" 


7 Be + 


e~ Li + v e 


" 8 B neutrinos" 


8 B -> f 


' Be* + e+ + v e 


"hep neutrinos" 


3 He + p 


^ 4 He + e+ + i/ e 



Table 1: Reactions from the Sun producing neutrinos. 



For low energy neutrino sources, the CC interaction may also be suppressed due 
to conversion of the nucleon at the lower vertex. For example, the CC interaction 
commonly called "inverse beta decay" (IBD), v e p — > e + n, which is crucial to reactor 
neutrino experiments, has a threshold of 1.084 MeV, driven by the mass difference 
between the proton and the neutron plus the mass of the positron. In the case of bound 
nuclei, the energy transferred in a CC interaction must overcome the binding energy 
difference between the incoming and outgoing nucleus as well as the mass suppression 
due to the charged lepton. This leads to nuclear-dependent thresholds for the CC 
interaction. For example: 



35 C1(75.8%) 


^35 


Ar : 


5.967 MeV; 


37 C1(24.2%) 


^37 


Ar : 


0.813 MeV; 


69 Ga(60.1%) - 


^69 


Ge : 


2.227 MeV; 


71 Ga(39.9%) - 


_^71 


Ge : 


0.232 MeV. 



are the thresholds for isotopes which have been used as targets in past solar neutrino 
(u e ) detectors. 14-16 

In discussing neutrino scattering at higher energies, several kinematic quantities 
are used to describe events. The squared center of mass energy is represented by the 
Mandelstam variable, s. The energy transferred by the boson is v, and y = vjE v is the 
fractional energy transfer, or "inelasticity." The distribution of events as a function of 
y depends on the helicity. For neutrino scattering from quarks, the y-dependence is 
flat, but for antineutrinos, the differential cross section is peaked at low y. The variable 
Q 2 is the negative squared four-momentum transfer. Deep inelastic scattering begins 
to occur at Q 2 ~ 1 GeV 2 . If x is the fractional momentum carried by a struck quark 
in a deep inelastic scatter, then x = Q 2 /2Mis, where M is the target mass. Elastic 
and quasieleastic scattering occur at x = 1, hence Q 2 = 2Mu ps sxy, valid for large s. 

1.3 Sources of Neutrinos 

With such a small interaction probability, it is clear that intense neutrino sources 
are needed to have high statistics in a neutrino experiment. The primary sources 
of neutrino for interactions observed on Earth are the Sun, cosmic-ray interactions, 
reactors, and accelerator beams. 

At present, there are two intense sources in the few MeV range that allow for low 
energy neutrino interaction studies. First, the interactions in the Sun produce a pure 
v e flux, as listed in tab. 1. The energy distribution of neutrinos produced by these 
reactions is shown in Fig. 4. The sensitivity of various solar neutrino experiments, 
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Neutrino Energy (MeV) 



Figure 4: The flux predicted by the Standard Solar Model. 13 The sensitivity of past solar 
neutrino detectors varies due to CC threshold in the target material. 14-16 The thresholds for 
various experiments is shown at the top of the plot. 

due to the CC threshold, is shown at the top of the figure. There is no observable 
antineutrino content. The best limit on the solar neutrino D e /v e ratio for E u > 8.3 
MeV is 2.8 x 10~ 4 at 90% CL. 17 The second source is from reactors. In contrast to the 
Sun, reactors produce a nearly pure v e flux. The energy peaks from ~ 3 to 7 MeV. 
Neutrinos from j3 decay of accelerated isotopes could, in principle, represent a third 
intense source of neutrinos in the MeV range (or higher), once the technical issues 
involved in designing such an accelerator are overcome. Such a "beta beam" would 
produce a very pure v e or v e beam, depending on the accelerated isotope. 18 

At present, higher energy experiments use neutrinos produced at accelerators and in 
the atmosphere. In both cases, neutrinos are dominantly produced via meson decays. 
In the atmospheric case, cosmic rays hit atmospheric nuclei producing a shower of 
mesons which may decay to neutrinos along their path through the atmosphere to 
Earth. In a conventional neutrino beam, protons impinge on a target, usually beryllium 
or carbon, producing secondary mesons. In many experiments, the charged mesons are 
focussed (bent) toward the direction of the experiment with a magnetic device called 
a horn. These devices are sign-selecting - they will focus one charge-sign and defocus 
the other - and so produce beams which are dominantly neutrinos or antineutrinos 
depending on the sign-selection. The beamline will have a long secondary meson decay 
region, which may be air or vaccuum. This is followed by a beam dump and an 
extended region of dirt or shielding to remove all particles except neutrinos. There 
are excellent reviews of methods of making accelerator-produced neutrino beams. 19 
Tab. 2 summarizes the common sources of neutrino production in the atmosphere and 
conventional accelerator based beams. 

Many atmospheric and accelerator-based neutrino experiments are designed to 
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2-body pion decay ir + — > n~ — > 

2-body kaon decay K + — > /Lt + z/ M , i^ - — > /U~z/ M 

muon decay /i + — ► e+z/^z/e, //~ — > e~v il p e 

K e3 decay lf + — > 7r°e + z/ e , X - — > 7T°e~z/ e , K° — > 7r~e + z/ e , iC° — > 7r + e~z/ e 

Table 2: Common sources of neutrinos in atmospheric and accelerator experiments. 

study 100 MeV to 10 GeV neutrinos. The atmospheric neutrino flux drops as a power- 
law with energy, and the 1 to 10 GeV range dominates the event rate. Accelerator 
beams can be tuned to a specific energy range and, using present facilities, can extend 
to as high as 500 GeV. From the viewpoint of sheer statistics, one should use the high- 
est energy neutrino beam which is practical for the physics to be addressed, since the 
cross section rises linearly with energy. However, lower neutrino energy beams, from 
~ 1 to 10 GeV, are typically used for oscillation experiments. In these experiments, 
having a cleanly identified lepton in a low multiplicity event trumps sheer rate, and so 
~ 1 GeV beams are selected to assure that CCQE and single pion events dominate the 
interactions. 

Both atmospheric and accelerator based neutrino sources are dominantly z^-flavor. 
The main source of these neutrinos is pion decay. To understand why pions preferen- 
tially decay to produce Un rather than v e , consider the case of pion decay to a lepton 
and an antineutrino: ir~ — ► £~i>£. The pion has spin zero and so the spins of the 
outgoing leptons from the decay must be opposite from angular momentum conserva- 
tion. In the center of mass of the pion, this implies that both the antineutrino and 
the charged lepton have spin projected along the direction of motion ( "right" or "pos- 
itive" helicity). However, this is a weak decay, where the W only couples to the RH 
antineutrino and the LH component of the charged particle. The amplitude for the 
LH component to have right-helicity is proportional to m/E. Thus it is very small for 
an electron compared to the muon, producing a significant suppression for decays to 
electrons. Calculating the expected branching ratios: 

- r ( 7r± ^ e±1/e ) (7\ 

^theory - r ^±_^± u ^ {< ) 

= {^)\ z<~S ) 2 ( 8 ) 

= 1.23 x 10~ 4 ; (9) 

This compares well to the data: 20 R exp = (1.230 ± 0.004) x 10" 4 . 

The above discussion assumed the neutrino was massless. If the neutrino is massive, 
then it too can be produced with wrong helicity with an amplitude proportional to 
m u /E and thus a probability proportional to (m u /E) 2 . As discussed in sec. 2.2, below, 
neutrino mass is limited to be very small (~ eV) and thus the rate of wrong- helicity 
neutrino production is too low a level for any chance of observation in the near future. 

Depending on the energy, there may also be significant neutrino production from 
kaon decays. The charged kaon preferentially decays to the for the same reason as 
the charged pion. However, for equal energy mesons, the kinematic limit for a neutrino 
from K + decay is much higher than for tt + decay: E™ ax ' K = 0.98£x compared to 
£jmax,n _ q Thus the neutrinos from kaon decays can be isolated by studying 

the highest energy component of a beam. Fig. 5 shows the contributions of pion and 
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Figure 5: The contributions from pion and kaon production to the total predicted flux in 
the MiniBooNE experiment. The spikes at low energy in the K-produced fluxes are due to 
decays of stopped kaons in the beam dump. 10 
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Figure 6: The variation in the atmospheric neutrino flavor content as a function of energy for 
two locations, Japan (solid line, open circles) and Minnesota (dashed line, closed squares). 
The points are from a full 3-dimensional monte carlo of the flux, while the histograms are 
from a simpler model. 21 

kaon decays to the flux in the MiniBooNE experiment, which uses an 8 GeV primary 
proton beam. 

Electron neutrino flavors are produced in these beams through K — > irv e e (called 
"Ke3" ) decay and through the decay of the muons which were produced in the pion 
decay. These are three-body decays which avoid substantial helicity suppression. He- 
licity does, however, affect the energy spectrum of the outgoing decay products. In an 
accelerator-based experiment, the level of electron-flavor content can be regulated, at 
some level, by the choice of primary beam energy and the length of the decay region. 
A low primary beam energy will suppress kaon production because of the relatively 
high mass of this meson (494 MeV) . A short decay pipe will suppress v e from fi decay, 
which tends to occur downstream, because it is produced in a multi-step decay chain 
(ir — > n — ► v e ). Both of these methods of suppressing v e production also lead to a 
reduction in the production rate, so an experimenter must balance competing goals 
in the beam design. In the case of atmospheric neutrinos, the ratio is roughly 2:1 
for Vp.Ve, though the fraction v e 's changes with energy (see Fig. 6). The atmospheric 
flux depends on the location of the detector because charged particle are bent by the 
Earth's magnetic field. The variation between fluxes at the Kamioka mine in Japan 
and Soudan mine in Minnesota are shown in Fig. 6. 

As we move to a precision era in neutrino physics, precise "first-principles" pre- 
dictions of the flux are becoming very important. For conventional accelerator-based 
neutrino beams and for the atmospheric flux, this requires well-measured cross-sections 
for production of secondary pions and kaons. This has motivated a range of secondary 
production experiments. The kinematic coverage is shown on Fig. 7. 

The future of high intensity z/„ and v e beams is likely to lie in beams produced from 
muon decay. Because of the potential for very high intensity, these beams are called 
"Neutrino Factories." The concept is very attractive because it produces beams which 
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Figure 7: The kinematic range covered by recent experiments measuring secondary pion and 
kaon production. 24 



are very pure Un and l> e from fi~ and vice versa from Each flavor has no "wrong 
sign" (antineutrino-in-neutrino-beam or neutrino-in-antineutrino beam) background. 
However, neutrino factory designs 18 necessarily produce high energy neutrinos, since 
the muons must be accelerated to high energies in order to live long enough to be 
captured and circulated in an accelerator. The Neutrino Factory is seen as a promising 
first machine for testing ideas for a muon collider, 18 and thus has attracted interest 
beyond the neutrino community. 

A beam enriched in v T can be produced by impinging very high energy protons on 
a target to produce D s -mesons which are sufficiently massive to decay to r v T . The 
r lepton is very massive, at 1.8 GeV, compared to the muon, at 106 MeV, and thus 
helicity considerations for the D s decay strongly favor the r v T mode compared to /i ^, 
by a ratio of about 10:1. The r then subsequently decays, also producing v T s. 

Unfortunately, because of the short lifetime, it is not possible to separate D s mesons 
from the other mesons prolifically produced by the primary interaction. As a result, 
the beam is dominated by the u^s produced by decays of other mesons. To reduce the 
production of v^, experiments use a "beam dump" design where protons hit a very 
thick target where pions can be absorbed before decaying. The only enriched-^ r beam 
created to date was developed by DoNuT. 6 They used an 800 GeV proton on a beam 
dump, to produce a ratio of v e :v ii :v T of about 6:9:1. 
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1.4 Typical Neutrino Detectors 



Because neutrinos interact so weakly, the options for detectors are limited to designs 
which can be constructed on a massive scale. There are several general styles in 
use today: unsegmented scintillator detectors, unsegmented Cerenkov detectors, seg- 
mented scintillator-and-iron calorimeters, and segmented scinitillator trackers. The 
most promising future technology is the noble-element based detector, which is effec- 
tively an electronic bubble chamber. Liquid argon detectors are likely to be the first 
large-scale working example of such technology. There are a few variations on these 
five themes, which are considered in later sections in the context of the measurement. 

Unsegmented scintillator detectors are typically used for low energy antineutrino 
experiments. Recent examples include Chooz, 25 KamLAND 26 and LSND. 27 These 
consist of large tanks of liquid scintillator surrounded by phototubes. Usually the scin- 
tillator is oil based, hence the target material is CH2 and its associated electrons. Often 
the tubes are in an pure oil buffer. This reduce backgrounds from radiation emitted 
from the glass which would excite scintillator. The free protons in the oil provide a 
target for the interaction, v e p — > e + n, which is the key for reactor experiments. The 
reaction threshold for this interaction is 1.806 MeV due to the mass differences between 
the proton and neutron and the mass of the positron. The scintillation light from the 
e + , as well as light from the Compton scattering of the 0.511 MeV annihilation pho- 
tons provide an initial ("prompt") signal. This is followed by n capture to produce 
deuterium and a 2.2 MeV. This sequence - positron followed by neutron capture - pro- 
vides a clean signal for the interaction. Doping the liquid scintillator with gadolinium 
substantially increases the neutron capture cross section as well as the visible energy 
produced in the form of gammas upon neutron capture. 

Unsegmented scintillator detectors are now being introduced for low energy so- 
lar neutrino measurements at Borexino, 28 KamLAND 29 and SNO+. 30 These provide 
energy information on an event-by-event basis, unlike most past solar neutrino exper- 
iments, such as Homestake, 14 SAGE 15 and GallEx, 16 which intergrated over time and 
energy. However, these are very difficult experiments to perform because a neutron is 
not produced and so the scattering does not produce a two-fold coincidence, but only 
a prompt flash of light. 

Environmental backgrounds are by far the most important issue in low energy 
experiments. These fall into two categories: naturally occurring radioactivity and 
muon-induced backgrounds. To get a sense for what is expected, Fig. 8 shows the 
visible energy distribution of singles events from the KamLAND experiment with the 
sources of environmental background identified. The naturally occurring radioactive 
contaminants mainly populate the low energy range of Fig. 8, with isotopes from the 
U and Th chain extending to the highest energies. These isotopes must be kept un- 
der control by maintaining very high standards of cleanliness. The second source of 
environmental background, the /3-decays of isotopes produced by cosmic ray muons. 
These dominate the background for E v iu e > 4 MeV (see Fig. 8). These can only be 
eliminated by shielding the detector from cosmic rays. As a result, we must build deep 
underground laboratories with many thousands of meters- water-equivalent ( "mwe" ) of 
rock shielding. 

In these scintillator detectors, the CC interaction with the carbon in the oil (which 
produces either nitrogen or boron depending on whether the scatterer is a neutrino 
or antineutrino) has a significantly higher energy threshold than scattering from free 
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Figure 8: Energy distribution and sources of singles events in KamLAND as a function of 
visible energy. 31 



Figure 9: An example of a muon ring (left) and electron ring (right) in the Super K Cerenkov 
detector. 33 

protons. v e + C — ► e~ + N has a threshold energy of 13.369 MeV, which arises from 
the carbon- nitrogen mass difference (plus the mass of the electron). In the case of both 
reactor and solar neutrinos, the flux cuts off below this energy threshold. 

Existing unsegmented Cerenkov detectors include MiniBooNE, 32 Super K, 33 and 
AMANDA. 34 These detectors make use of a target which is a large volume of a clear 
medium (undoped oil, water and ice, respectively) surrounded by or interspersed with 
phototubes. Undoped oil has the advantages of a larger refractive index, leading to 
larger Cerenkov opening angle, and of not requiring a purification system to remove 
living organisms. Water is the only affordable medium once a detector is larger than a 
few ktons. For ultra-high energy neutrino experiments, a vast natural target is needed. 
Sea water 35 and ice 34 have been used. Ice is, to date, more successful because it does 
not suffer from backgrounds from bioluminescence. 

In most cases of these detectors, the tubes surround the medium and the projected 
image of the Cerenkov ring is used for particle identification. To understand how 
this works, first consider the case of a perfect, short track. This will project a ring 
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with a sharp inner and outer edge onto the phototubes. Next consider an electron 
produced in a v e CC quasielastic interaction. Because the electron is low mass, it will 
multiple scatter and easily bremsstrahlung, smearing the light projected on the tubes 
and producing a "fuzzy" ring. A muon produced by a CC quasileastic f M interaction 
is heavier and thus will produce a sharper outer edge to the ring. For the same visible 
energy, the track will also extend farther, filling the interior of the ring, and perhaps 
exit the tank. Fig. 9 compares an electron and muon ring observed in the Super K 
detector. If the muon stops within the tank and subsequently decays, the resulting 
"michel electron" provides an added tag for particle identification. In the case of the 
18% will capture in water, and thus have no michel electron tag, while only 8% 
will capture in oil. 

Scintillator and iron calorimeters provide affordable detection for interactions in 
the range of ~1 GeV and higher. Recent examples include the MINOS 36 and NuTeV 37 
experiments. In these detectors, the iron provides the target, while the scintillator 
provides information on energy deposition per unit length. This allows separation be- 
tween the hadronic shower, which occurs in both NC and CC events, and the minimum 
ionizing track of an outgoing muon, which occurs in CC events. Transverse informa- 
tion can be obtained if segmented scintillator strips are used, or if drift chambers are 
interspersed. The light from scintillator strips is transported to tubes by mirrored 
wave- length-shifting fibers. Transverse information improves separation of electromag- 
netic and hadronic showers. The iron can be magnetized to allow separation of neutrino 
and antineutrino events based on the charge of the outgoing lepton. 

In all three of the above detector designs, it is difficult to reconstruct multi-particle 
events. Tracking is not an option for an unsegmented scintillator detector. Cerenkov 
detectors can typically resolve two tracks per event. Segmented calorimeters reduce 
multiple hadrons to a shower, obscuring any track-by-track information other than 
from muons. 

To address the problem of track reconstruction in low energy (< 1 GeV), low 
multiplicity events, there has been a move toward all-scintillator tracking detectors. 
This began with the SciBar detector in K2K. 38 This detector used scintillator strips, 
as in MINOS, but without interspersing iron. As a result, low energy (few MeV) 
tracks were clearly observable and quasielastic and single pion events could be fully 
reconstructed. SciBar has since been incorporated into the SciBooNE experiment at 
Fermilab. 11 The CCQE event in SciBooNE, shown in Fig. 10, makes clear the benefits 
of fine segmentation. The position of the vertex and the short track from the proton 
are well-resolved in the SciBar detector (green region). The technology has been taken 
further by the MINERvA experiment, which has attained 2 mm resolution with their 
prototype. 39 Scibar and MINERvA are relatively small (few ton) detectors. The 
first very large scale application of this technology will be NOvA, which is a future 
15 kton detector. 40 This detector will use PVC tubes filled with liquid scintillator, 
which is more cost-effective than extruded scintillator strips for very large detectors. 
Their design also loops the wave-length shifting fiber, so that there are, effectively, two 
perfectly mirrored fibers are in each cell. This elegant solution increases the collected 
light by a factor of four, which is necessary for ~ 15 m strips. 

The most promising new technology for high resolution track reconstruction in 
neutrino physics is the liquid argon TPC. A TPC, or time projection chamber, uses 
drift chambers to track in the x and y views and drift time to determine the z view. 
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Figure 10: A CCQE(z/^ + n — > /x + p) event observed in the SciBooNE detector. The long, 
minimum- ionizing red track is identified as the muon, the short, heavily- ionizing red track 
is identifed as the proton. 11 



Liquid argon (LAr), which provides the massive target for the neutrino interaction, also 
scintillates, providing the start for the drift-time measurement. A key point for future 
neutrino experiments is the high efficiency for identifying electron showers (expected to 
be 80-90%) with a rejection factor of 70 for NC tt° events. In particular, these detectors 
can differentiate between converted photons and electrons through the dE/dx in the 
first few centimeters of the track. Typical energy resolution for an electromagnetic 
shower is 3%/V~E. 

There is a great deal of activity on development of LAr detectors. Data have been 
taken successfully on a 50 liter LArTPC prototype in the NOMAD neutrino beam at 
CERN, resulting in reconstruction of ~100 CC quasielastic events. 41 Also, recently, a 
600 ton Icarus module has been commissioned at Gran Sasso. 42 A 0.8 ton LAr test 
detector will begin taking data at Fermilab in January, 2008. 43 As discussed in sec. 3.2, 
the microBooNE experiment is a proposed 100 ton detector which would take data in 
20 1 0. 44 In principle, these detectors can be scaled up to tens of ktons, as is discussed 
in the "Ash River Proposal" . 45 



2 Neutrinos As We Know Them Now 

The recent discovery of neutrino oscillations requires that we reconsider the Standard 
Model Lagrangian of sec. 1.1. It must now incorporate, preferably in a motivated 
fashion, both neutrino mass and neutrino mixing. This represents both a challenge 
and an opportunity for the theory, which I will discuss in the following section. This 
section concentrates on the experimental discovery. It is interesting to note that while 
neither neutrino mass nor mixing were "needed" in the Standard Model theory, both 
are required for the discovery of neutrino oscillations. The probability for neutrino 
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Figure 11: Example of neutrino oscillations as a function of distance from the source, L. The 
wavelength depends upon the experimental parameters L and E (neutrino energy) and the 
fundamental parameter Am 2 . The amplitude of the oscillation is constrained by the mixing 
term, sin 2 29. 

oscillations will be zero unless both effects are present. 

The outcome of the observation of neutrino oscillations is typically summarized 
by the statement that "neutrinos have mass." To be clear: we still have no direct 
measurement of neutrino mass. At this point, we have clear evidence of mass differences 
between neutrinos from the observation of neutrino oscillations. A mass difference 
between two neutrinos necessarily implies that at least one of the neutrinos has non- 
zero mass. All experimental evidence indicates that the actual values of the neutrino 
masses are tiny in comparison to the masses of the charged fermions. At the end of 
this section, attempts at direct measurement of neutrino mass are described. 

2.1 Neutrino Oscillations 

Recent results on neutrino oscillations provide indisputable evidence that there is a 
spectrum of masses for neutrinos. In this section, I describe the formalism for neutrino 
oscillations, and then review the experimental results which have now been confirmed 
at the 5a level. This is covered briefly because these results are well known and covered 
extensively elsewhere. 2 

2.1.1 The Basic Formalism 

Neutrino oscillations requires that neutrinos have mass, that the difference between the 
masses be small, and that the mass eigenstates be different from the weak interaction 
eigenstates. In this case, the weak eigenstates can be written as mixtures of the mass 
eigenstates. For example, in a simple 2-neutrino model: 

v e = cos 9 v\ + sin 9 vi 
Vfj, = — sin 9 v\ + cos 9 v<i 
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where 6 is the "mixing angle." In this case, a pure flavor (weak) eigenstate born through 
a weak decay can oscillate into another flavor as the state propagates in space. This 
oscillation is due to the fact that each of the mass eigenstate components propagates 
with different frequencies if the masses are different, Am 2 = |m| — m\\ > 0. In such a 
two-component model, the oscillation probability for — > v e oscillations is then given 
by: 

9 9 A-27 Am 2 (eV 2 ) L(km)\ , 
Prob [u, - v e ) = sin 2 20 sin 2 ( E ^ Q J ) j , (10) 

where L is the distance from the source, and E is the neutrino energy. As shown in 
Fig. 11, the oscillation wavelength will depend upon L, E, and Am 2 . The amplitude 
will depend upon sin 2 29. 

Neutrino oscillations only occur if the two mass states involved have sufficiently 
small Am 2 that the neutrino flavor is produced in a superposition of two mass states. 
If the mass splitting is sufficiently large, a given neutrino flavor would be produced in 
one or the other of the two mass eigenstates and interference (i.e., oscillations) would 
not occur. 

Most neutrino oscillation analyses consider only two-generation mixing scenarios, 
but the more general case includes oscillations among all three neutrino species. This 
can be expressed as: 







( u el 


u e2 


u e3 




H 






u& 


\"t J 




V U T 1 


U r2 


u T3 







II 









This formalism is analogous to the quark sector, where strong and weak eigenstates are 
not identical and the resultant mixing is described conventionally by a unitary mixing 
matrix. The oscillation probability is then: 

Prob (v a -> vp) = S a p - 

j>i \ J 

where Am 2 ^ = m 2 — mf, a and are flavor-state indices (e,/x, r) and i and j are 
mass-state indices (1,2,3). 

For three neutrino mass states, there are three different Am 2 parameters, although 
only two are independent since the two small Am 2 parameters must sum to the largest. 
The neutrino mass states, vi, v<i and v 3 are defined such that the difference between v\ 
and vi always represents the smallest splitting. However, the mass of v 3 relative to v\ 
and V2 is arbitrary and so the sign of the Am 2 parameters which include the third mass 
state may be positive or negative. That is, if v 3 > ^1,^2, then Am| 3 will be positive, 
but if v\,vi > 1/3, then Am| 3 will be negative. The former is called a "normal mass 
hierarchy" and the latter is the "inverted mass hierarchy." At this point, the sign is 
irrelevant because Am 2 appears in a term which is squared. However, in sec. 3.1, this 
point will become important. 

The mixing matrix above can be described in terms of three mixing angles, #12, #13 
and 6*23: 

C12C13 S12C13 S13 \ 

U = I -S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 I , (12) 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13/ 
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where Cjj = cos dij and Sij = sin with iand j referring to the mass states. In fits 
to the oscillation parameters, people variously quote the results in terms of the matrix 
element of U, sin-squared of the given angle, sin-squared of twice the angle and a 
variety of other forms, all of which are related. Using the 13 case as an example, the 
quoted parameters are related by: 



Thus, in total, there are five free parameters in the simplest three-neutrino oscilla- 
tion model, which can be taken to be Am 2 2 , Ara^, #12, #13 and 623- 

Although in general there will be mixing among all three flavors of neutrinos, two- 
generation mixing is often assumed for simplicity. If the mass scales are quite different 
(1713 » 771,2 >> Tni, for example), then the oscillation phenomena tend to decou- 
ple and the two-generation mixing model is a good approximation in limited regions. 
In this case, each transition can be described by a two-generation mixing equation. 
However, it is possible that experimental results interpreted within the two-generation 
mixing formalism may indicate very different Am 2 scales with quite different appar- 
ent strengths for the same oscillation. This is because, as is evident from equation 11, 
multiple terms involving different mixing strengths and Am 2 values contribute to the 
transition probability for u a — > v$. 

2.1.2 Matter Effects 

The probability for neutrino oscillations is modified in the presence of matter. This 
is true in any material, however the idea was first explored for neutrino oscillations in 
the Sun, by Mikheyev, Smirnov and Wolfenstein. Therefore, matter effects are often 
called "MSW" effects. 46 In general, matter effects arise in neutrino-electron scattering. 
The electron neutrino flavor experiences both CC and NC elastic forward-scattering 
with electrons. However, the v„ and v T experience only NC forward-scattering, because 
creation of the fi or r is kinematically forbidden or suppressed (e.g. u^ + e~ — > v e +[i~). 
This difference produces the matter effect. 

For neutrinos propagating through a constant density of electrons, if V e is the elastic 
forward scattering potential for the v e component, and V ther is the potential for the 
other neutrino flavors, then the additional scattering potential is 



where Gp is the Fermi constant and n e is the electron density. This potential mod- 
ifies the Hamiltonian, so that, if Hq is the vacuum Hamiltonian, then in matter the 
Hamiltonian is Hq + V. This means that the eigenstates are modified from those of a 
vacuum, v\ and 1/2, to become v\ m and V2m- Effectively, the neutrino mass spectrum is 
not the same as in vacuum. The solutions to the Hamiltonian are also modified. From 
this, one can see that the presence of electrons may substantially change the oscillatory 
behavior of neutrinos. 

The simplest outcome is that matter induces a shift in the mass state, which is a 
combination of flavor eigenstates, propagates through the material. This leads to a 
change in the oscillation probability: 



U 2 3 ps sin 2 6i3 ps -sin 2 29i3. 



(13) 



V = V e - V ther = V2G F n e , 



(14) 



Prob (v e — > v t 



7O = (sin 2 20/W 2 ) sin 2 (l.27WAm 2 



LjE) 



(15) 
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where W 2 = sin 2 26 + {\/2G pn e (2E / Am 2 ) — cos26>) 2 (Note that in a vacuum, where 
n e = 0, this reduces to equation 10.) From this, one can see that if a neutrino, 
passing through matter, encounters an optimal density of electrons, a "resonance," or 
large enhancement of the oscillation probability, can occur. The Sun has a wide range 
of electron densities and thus is a prime candidate for causing matter effects. Also, 
neutrinos traveling through the Earth's core, which has a high electron density, might 
experience matter effects. This will produce a "day-night effect," or siderial variation, 
for neutrinos from the Sun. 

For situations like the Sun, with very high electron densities which vary with the 
position of the neutrino (and hence the time which the neutrino has lived), the situation 
is complex. If the electron density is high and the density variation occurs slowly, or 
adiabatically, then transition (not oscillation!) between flavors in the mass state can 
occur as the neutrino propagates. Thus it is possible for neutrinos to be produced 
in the core of the Sun in a given mass and flavor state, and slowly evolve in flavor 
content until the neutrino exits the Sun, still in the same mass state. In other words, 
in the Sun, a v e produced in a mass eigenstate V2 m (r), which depends on the local 
electron density at radius r, propagates as a vim until it reaches the r = R so iar, where 
V2m(Rsoiar) = ^2- This peculiar effect is called the Large Mixing Angle MSW solution. 

2.1.3 Designing an Oscillation Experiment 

From equation 10, one can see that three important issues confront the designer of the 
ideal neutrino experiment. First, if one is searching for oscillations in the very small 
Am 2 region, then large L/E must be chosen in order to enhance the sin 2 (1.27Am 2 L/£') 
term. However if L/E is too large in comparison to Am 2 , then oscillations occur 
rapidly. Because experiments have finite resolution on L and E, and a spread in beam 
energies, the sin 2 (1.27Am 2 L/£') averages to 1/2 when Am 2 S> L/E and one loses 
sensitivity to Am 2 . Finally, because the probability is directly proportional to sin 2 26, 
if the mixing angle is small, then high statistics are required to observe an oscillation 
signal. 

There are two types of oscillation searches: "disappearance" and "appearance." 
To be simplistic, consider a pure source of neutrinos of type x. In a disappearance 
experiment, one looks for a deficit in the expected flux of v x . This requires accurate 
knowledge of the flux, which is often difficult to predict from first principles. Therefore, 
most modern disappearance experiments employ a near-far detector design. The near 
detector measures the flux prior to oscillation (the design goal is to effectively locate 
it at L = in Fig. 11). This is then used to predict the unoscillated event rate in the 
far detector. A deficit compared to prediction indicates disappearance. Appearance 
experiments search for v a — > v$ by directly observing interactions of neutrinos of type 
(3. The case for oscillations is most persuasive if the deficit or excess has the (L/E) 
dependence predicted by the neutrino oscillation formula (equation 10). 

The "sensitivity" of an experiment is defined as the average expected limit if the 
experiment were performed many times with no true signal (only background) . Let us 
consider the sensitivity for a hypothetical perfect (no-systematic error) disappearance 
neutrino oscillation experiment with N events. A typical choice of confidence level is 
90%, so in this case, the limiting probability, assuming there is no signal, is 

P = aVN/N. (16) 
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Figure 12: An illustration of the sensitivity of an imaginary oscillation experiment. The 
region of sensitivity for an experiment depends on the oscillation probability, P, where one 
can set a limit at some confidence level. Most experiments use 90% CL. The boundaries 
depend on P, L and E. 

There are two possible choices of a associated with a 90% CL sensitivity, depending 
on the underlying philosophy. If one assumes there is no signal in the data, then one 
quotes the sensitivity based on 90% of a single-sided Gaussian, which is a = 1.28. 
If the philosophy is that there is a signal which is too small to measure, then one 
quotes the sensitivity using a = 1.64, which is appropriate for a double-sided Gaussian. 
Historically, a = 1.28 was used in most publications. Physicists engage in arguments 
as to which is most correct, but what is most important from a practical point of view 
is for the reader to understand what was used. The reader can always scale between 
1.28 and 1.64 depending on personal opinion. 

There is only one measurement, P, and there are two unknowns, Am 2 and sin 2 29; 
so this translates to a region of sensitivity within Am 2 - sin 2 29 space. This is typically 
indicated by a solid line, with the allowed region on the right on a plot (see illustration 
in Fig 12). For the perfect (no-systematic error) experiment, the high Am 2 limit on 
sin 2 29 is driven by the statistics. On the other hand, the L and E of the experiment 
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drive the low Am 2 limit, which depends on the fourth root of the statistics. If our 
perfect experiment had seen a signal, the indications of neutrino oscillations would 
appear as "allowed regions," or shaded areas on plots of Am 2 vs. sin 2 26. 

This rule of thumb - that statistics drives the sin 2 2#-reach and L/E drives the Am 2 
reach - becomes more complicated when systematics are considered. The imperfections 
of a real experiment affect the limits which can be set. Systematic uncertainties in the 
efficiencies and backgrounds reduce the sensitivity of a given experiment. Background 
sources introduce multiple flavors of neutrinos in the beam. Misidentification of the 
interacting neutrino flavor in the detector can mimic oscillation signatures. In addition, 
systematic uncertainties in the relative acceptance versus distance and energy need to 
be understood and included in the analysis of the data. 

For a real experiment, with both statistical and systematic errors, finding the sensi- 
tivity and final limit or allowed region requires a fit to the data. The data are compared 
to the expectation for oscillation across the range of oscillation parameters, and the set 
of parameters where the agreement is good to 90% CL are chosen. Historically, there 
are three main approaches which have been used in fits. The first method is the "single 
sided raster scan." In this case one chooses a Am 2 value and scans through the sin 2 26- 
space to find the 90% CL limit. The second method, the "global scan," explores Am 2 - 
and sin 2 26>-space simultaneously. Thus there are two parameters to fit and two degrees 
of freedom. The third method is the frequentist, or "Feldman-Cousins" approach, 47 
in which one simulates "fake-experiments" for each Am 2 and sin 2 26 point, and deter- 
mines the limit where, in 90% of the cases, no signal is observed. Each method has 
pros and cons and the choice is something of a matter of taste. As with the question of 
a single- or double-sided gaussian, what is important is to compare sensitivities, limits, 
and signals from like methods. 

It is possible for an experiment which does not observe a signal to set a limit which 
is better than the sensitivity. This occurs if the experiment observed a downward 
fluctuation in the background. In this limit is hard to interpret. The latest 

standard practice is to show the sensitivity and the limit on plots, and the readers can 
draw their own interpretation. 47 

2.1.4 Experimental Evidence for Oscillations 

Two separate allowed regions in Am 2 -and-sin 2 29 -space for neutrino oscillations have 
been observed at the > 5cr level. These are called the "Atmospheric Am 2 " and "Solar 
Am 2 " regions. The names are historical, as will be seen below. Many reviews have 
been written on these results (see, for example,, 2 , 48 and 49 ) and so here the results are 
briefly outlined. 

The highest Am 2 signal was first observed using neutrinos produced in the upper 
atmosphere. These atmospheric neutrinos are produced through collisions of cosmic 
rays with the atmosphere. The neutrinos are detected through their charged-current 
interactions in detectors on the Earth's surface. 

The first evidence for atmospheric neutrino oscillations came from the Kamioka 50 
and 1MB 51 experiments. This was followed by the convincing case presented by the 
Super K experiment. 53 These were single detector experiments observing atmospheric 
neutrino interactions as a function of zenith angle (see Fig. 13). Several striking features 
were observed. The first was that the flavor neutrinos showed clear evidence of 
disappearance while the v e flavor CC scatters were in good agreement with prediction. 
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Figure 13: Event rates observed in Super K as a function of zenith angle for two energy 
ranges. Candidate v e events are on the left, are on the right. The red line indicates the 
predicted rate. The green line is the best fit including oscillations. 52 



The second striking observation was that the apparent mixing was nearly maximal. 
In other words, the experiments were seeing a 50% reduction of the event rate 
compared to expectation. 

Complications in the analysis arise from the difficulty in understanding production 
of atmospheric neutrinos (affecting the understanding of E) and in the accurate re- 
construction of events as a function of zenith angle (affecting knowledge of L). The 
Am 2 extracted from the Kamoiokande data is an order of magnitude higher than that 
extracted from the Super K data, indicating a clear systematic effect. Thus, it was ab- 
solutely crucial for accelerator-based "long-baseline" neutrino experiments to confirm 
this result. In these experiments, the L is well defined by the distance from source to 
detector, and the E is well understood from a near detector measurement. 

The challenge for long-baseline experiments is that the L/E required to access the 
atmospheric signal is on the order of 1000 km/GeV. If the beam is relatively low energy, 
so that the easy-to-reconstruct CCQE interaction dominates the events, then L is on 
the order of 1000 km. This leads to two major technical challenges. First, because 
the Earth is a sphere, if the source and detector are to be located on (or near) the 
surface, the beam must be directed downward, into the Earth. Engineering a beamline 
at a steep angle requires overcoming substantial hurdles in tunneling. Second, the 
beam spreads as it travels outward from the source, resulting in low intensity at the 
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detector. Therefore, very high rates are needed. However, these challenges have now 
been overcome at three accelerator complexes: KEK, FNAL and CERN, and a new 
long-baseline beam from the JPARC facility will be available soon. Making use of these 
lines, initial confirmation of the atmospheric neutrino deficit came from the KEK-to- 
Kamiokande (K2K) long baseline experiment. 54 This has since been followed up by 
the MINOS experiment to high precision. 55 ' 56 

In the atmospheric data, the v e CC signal is in agreement with expectation, and 
in the long-baseline experiments, no v e excess has been observed. Therefore, one 
cannot interpret this oscillation signal as — > v e . This leaves only — ► v T as an 
explanation for the deficit in a three-neutrino model. Observation of v T CC interactions 
is experimentally difficult in these experiments for a number of reasons. First, the L/E 
of the signal is such that for lengths available to present experiments, the energy of 
the beam must be low (< 10 GeV). As discussed in sec. 1.2, because the r mass is 1.8 
GeV, there is substantial mass suppression for r production at low energies, so the CC 
event rate is low. Second, the r decays quickly, leaving behind a complicated event 
structure which can be easily confused with and v e low multiplicity interactions 
in calorimeter or Cerenkov detectors. The difficulty of identifying v T events even in 
a specialized emulsion-based detector with a high energy neutrino beam, was made 
clear by the DoNuT experiment, 6 which provided the first, and so far only, direct 
observation of CC v T interactions. Thus, while some studies claim observation of u T 
CC interactions in SuperK, 57 these results are not very convincing to this author. 
Fortunately, a specialized experiment called OPERA, 58 which is an emulsion-based 
long-baseline detector, is presently taking data. The average energy of the CNGS beam 
used by this experiment is 17 GeV, sufficiently high to produce v T CC events. This 
experiment is expected to observe ~ 15 events in 5 years of running if the atmospheric 
neutrino deficit is due to f« — ► v T with Am 2 = 2.5 x 10 -3 eV 2 . 59 

The lower Am 2 signal is called the "Solar Neutrino Deficit," as it was first observed 
as a low rate of observed z/ e 's from the Sun. The first observation of this effect was a 
v e deficit observed using a CI target 14 by Ray Davis and collaborators at Homestake, 
using v e + CI — ► e + Ar. Only about 1/3 of the total expected neutrino event rate 
was observed. By 1999, four additional experiments had confirmed these observations. 
The GALLEX 16 and SAGE 15 experiments confirmed a deficit for CC electron neutrino 
interactions in a Ga target producing Ge. The Super Kamiokande experiment observed 
a deficit for v e + e — > v e + e reactions in water. 22 The deficit is shown on Fig. 14, 
indicated by the blue points. This plot shows the ratio to the Standard Solar Model 
prediction, which is indicated by the solid line at unity. 

A few aspects of the initial solar neutrino deficit studies should be noted. First, the 
three types of experiments, chlorine-based, gallium-based, and water-based, measured 
different levels of deficit. Given that each type of nucleus has a different low energy 
threshold for observation of CC events, as previously discussed, one can interpret the 
varying levels of deficit as an energy dependent effect. Second, all of the above exper- 
iments rely upon the CC interaction. The energy of neutrinos from the Sun is so low, 
that should Vn or v T be produced through oscillations, the CC interaction could not 
occur. This is because of the relatively high mass of the \i (106 MeV) and the r (1.8 
GeV). Thus all of these experiments can observe that u e s disappeared, but they cannot 
observe if the neutrinos reappear as one of the other flavors. This makes a decisive 
statement that the effect is due to neutrino oscillations problematic. 
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Figure 14: Ratio of observed event rates in solar neutrino experiments compared to the 
Standard Solar Model. Experiments are plotted at the average energy of the detected signal, 
which varies due to detection threshold. Black error bars indicate Standard Solar Model 
error. 

For some time, people argued the apparent deficit was due to an incomplete picture 
of solar processes. The two important theoretical issues related to the solar neutrino 
fluxes were the fusion cross sections and the temperature of the solar interior. A 
comprehensive analysis of the available information on nuclear fusion cross sections 
important to solar processes has been compiled 60 and shows that the important cross 
sections are well-known. Results in helioseismology provided an important further 
test of the "Standard Solar Model" . 61 The Sun is a resonant cavity, with oscillation 
frequencies dependent upon P/p, the ratio of pressure to density. Helioseismological 
data confirmed the SSM prediction of U to better than 0.1%. 62 With the results 
of these studies, most physicists were convinced that the Standard Solar Model was 
substantially correct. The error bars on the black line at unity in Fig. 14 shows the 
side of the estimated systematic error on the Standard Solar Model. 

Interpreting the results as neutrino oscillations resulted in a complicated picture. 
The vacuum oscillation probability, calculated using equation 10, results in allowed 
regions of Am 2 which are very low (Am 2 ~ 10~ 10 eV 2 ). This is because the energy 
of the neutrinos is only a few MeV, and the Sun to the Earth pathlength is very long 
(~ 10 11 m) . On the other hand, the Sun has high electron content and density, so 
matter effects (sec. 2.1.2) could interfere with the picture, allowing higher true values 
of Am 2 . The MSW effect yielded two solutions in fits to the data. One was at mixing 
angles of ~ 10~ 3 . Until very recently, this was regarded as the most likely solution 
based on analogy with mixing in the quark sector. The other solution gave a very 
large, although not maximal, mixing angle. 

Two dramatic results of the early 2000 's demonstrated that the solar neutrino deficit 
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Figure 15: Events in KamLAND as a function of energy. The grey line indicates the expec- 
tation for no oscillation. 66 

was due to oscillations with the MSW effect and with large mixing angle. The first 
result was from the SNO experiment. 63 SNO used a D2O target which allowed for mea- 
surement of both CC v e interactions as well as u+d — ► u+n+p. In the former measure- 
ment, SNO sees a deficit consistent with the other measurements within an oscillation 
interpretation, and which yields a v e flux of (1.76±0.05(stat) ±0.09(sys)) x 10 6 /cm 2 s. 
The later measurement is an NC interaction, and thus is flavor-blind. It yields a to- 
tal NC flux of (5.09[}j!(stat)+o;4 6 (sys)) x 10 6 /cm 2 s 64 which can be compared with the 
theoretical prediction of (5.69 db 0.91) x 10 6 /cm 2 s. 65 In other words, SNO observed 
the expected total event rate, within errors. This implied that the v e s are oscillating 
to neutrinos which participate in the NC interaction, i/„s and/or v T s, with the total 
Vp + Vr flux equal to (3.41 ±0.45(stat)+°;^(sys)) x 10 6 /cm 2 s. 64 The results of two runs 
of the SNO experiment are shown by the red and green points on Fig. 14. The sec- 
ond result was from the KamLAND experiment. This was a reactor-based experiment 
located in Japan. Using many reactors which were hundreds of kilometers away, the 
KamLAND experiment was able to reach L/E ~ 1CT 6 m/MeV. This covered the MSW 
allowed-Am 2 solution. The statistics were on the order of hundreds of events, but this 
was enough to probe the large mixing-angle MSW solution. KamLAND expected 365 
events and observed 258 events, and thus had clear evidence for oscillations with large 
mixing, tan 2 9 = 0A0±° ]° and relatively high Am 2 , of 7.9±g;| x 1(T 5 eV 2 . 66 The 
energy distribution of the events observed in KamLAND is shown in Fig. 15. 

Based on the atmospheric and solar studies, there are two squared mass differences: 
Arnj. olar and Am\ tmos . The smaller is identified with the mass splitting between v\ and 
V2~. Am 2 2 = Am? solar . The atmospheric deficit measures a combination of Am^ and 
Amf 3 . However, since Am 2 3 = Am 2 2 + Am| 3 and Amf 2 is small, Am 2 3 m Am 23 ~ 

/ -^ m atmos- 
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Model(s) 


Refs. 


approximate sin 2 2# 13 


Minimal SO(10) 


68 


0.13 


Orbifold SO (10) 


69 


0.04 


SO(10) + Flavor symmetry 


70 


1.2 • 10~ 6 




71 


7.8 • 10" 4 




72-74 


0.01 .. 0.04 




75-77 


0.09 .. 0.18 


SO(10) + Texture 


78 


4- 10" 4 .. 0.01 




79 


0.04 


SU(2) L x SU(2) R x SU(4) C 


80 


0.09 


Flavor symmetries 


81-83 







84,85,94 


< 0.004 




87-89 


10~ 4 .. 0.02 




90-94 


0.04 .. 0.15 


Textures 


95 


4- 10~ 4 .. 0.01 




96—99 


0.03 .. 0.15 


3x2 see-saw 


100 


0.04 


Anarchy 


101 


> 0.04 


Renormalization group enhancement 


102 


0.03 .. 0.04 


M-Theory model 


103 


10" 4 



Table 3: Selected predictions for sin 2 26 



A recent global analysis of the data from the above experiments yields a consistent 
picture for three neutrino oscillations with five free parameters. The mass differences 
are: Amf 2 = (7.9 ± 0.3) x 10" 5 eV 2 and |Am 2 3 | = (2.5±g;|g) x 10" 5 eV 2 , where the 
absolute value is indicated in the second case because the sign (i.e. the mass hierarchy) 
is unknown. The two well-measured mixing angles are determined to be: sin 2 612 = 
0.30^q'q3 and sin 2 #23 = 0.507jo'o7- One mixing angle, 6*13 is yet to be measured, but a 
limit of sin 2 6*13 < 0.025 can be placed based on global fits. 

Based on the measurements, the mixing matrix of eq. 12, translates roughly into: 

/ 0.8 0.5 ? \ 
U = 0.4 0.6 0.7 . (17) 
\ 0.4 0.6 0.7 / 

This matrix, with its large off-diagonal components, looks very different from the quark- 
sector mixing matrix where the off-diagonal elements are all relatively small. In this 
matrix, the "odd element out" is U e 3 which is clearly substantially smaller than the 
others. At this point, there is no consensus on what this matrix may be telling us 
about the larger theory, but there is a sense that the value of 6*13 is an important 
clue. Theories which attempt to explain this matrix tend to fall into two classes - 
those where #13 is just below the present limit and those with very small values. As an 
illustration of this point, Tab. 3 shows order of magnitude predictions for a variety of 
theories. Thus a measurement of sin 2 2^13 which is greater than about 1%, or a limit 
at this level, can point the way to the larger theory. 
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The best method for measuring #13 is from reactor experiments which constrain 
this oscillation by searching for v e disappearance. The oscillation probability is given 
by: 

Preactor ^ sin 2 20 i3 sin 2 A + a 2 A 2 cos 4 0i3 sin 2 20i2, (18) 

with 

a = Ama/Ara^ (19) 
A = Am^L/(4£,). (20) 

Events are detected through the inverse beta decay (IBD) interaction. The CHOOZ 
experiment, 25 with a baseline of 1.1 km and typical neutrino event energies between 3 
and 5 MeV (E) = 3.5 MeV) has set the best reactor-based limit to date, of sin 2 20i 3 < 
0.27 at Am 2 = 2.5 x 10~ 3 eV 2 . This limit can be improved with a global fit, as quoted 
above. 

Significant improvement is expected from the upcoming round of reactor experi- 
ments results due to introducing a near-far detector design. The near detector mea- 
sures the unoscillated event rate, and the far detector is used to search for a deficit as 
a function of energy. The Double Chooz experiment, beginning in 2009, is expected 
to reach sin 2 20i 3 ~ 0.03. 105 This will be followed by the Daya Bay experiment which 
will reach ~0.01. 106 



2.2 Direct Measurements of Neutrino Mass 

For neutrinos, there are no mass measurements, only mass limits. Observations of 
neutrino oscillations are sensitive to the mass differences between neutrinos, not the 
actual mass of the neutrino. Therefore, they do not fall into the category of a "direct 
measurement". One can, however, use these oscillation results to estimate the required 
sensitivity for a direct mass measurement. The upper limit comes from assuming that 
one of the neutrino masses is exactly zero. Given that the largest Am 2 is ~ 3 x 10~ 3 
eV 2 , this implies there is a neutrino with mass V Am 2 ~ 0.05 eV. The mass of the 
neutrino can be directly measured from decay kinematics and from time of flight from 
supernovae. Neither method has reached the 0.05 eV range yet, although the next 
generation of decay-based experiments comes close. 

We know from neutrino oscillations that there is a very poor correspondence be- 
tween neutrino flavors and neutrino masses, i.e., the mixings are large. However, it is 
easiest to conduct the discussion of these limits in terms of specific flavors. Thus, what 
is actually being studied is an average mass associated with each flavor. For example, 
for the v e mass measured from (5 decay, which will be expanded upon below, what is 
actually probed is: 

m p = ^ t \U e ^m 2 . (21) 

The simplest method for measuring neutrino mass is applied to the v^. The mass 
is obtained from the 2-body decay-at-rest kinematics of it — > jiv^. One begins in the 
center of mass with the 4- vector relationship: p n = p^ + p v . Squaring and solving for 

neutrino mass gives: m 2 = m 2 + m 2 — ^/4m 2 r (|p M | 2 + m 2 ). From this, one can see that 
this technique requires accurate measurement of the muon momentum, p M , as well as 
the masses of the muon, m M and the pion, m n . In fact, the uncertainty on the mass 
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Experiment 


measured m 2 (eV 2 ) 


limit (eV), 95% C.L. 


Year 


Mainz 110 


-0.6± 2.2± 2.1 


2.2 


2004 


Troitsk 111 


-1.0 ± 3.0± 2.1 


2.5 


2000 


Mainz 112 


-3.7 ± 5.3 ± 2.1 


2.8 


2000 


LLNL 113 


- 130 ± 20± 15 


7.0 


1995 


CIAE 114 


- 31 ± 75± 48 


12.4 


1995 


Zurich 115 


-24 ± 48± 61 


11.7 


1992 


Tokyo INS 116 


- 65 ± 85± 65 


13.1 


1991 


Los Alamos 117 


- 147 ± 68± 41 


9.3 


1991 



Table 4: Overview of v e squared mass measurements. 



of the pion is what dominates the f M mass measurement. As a result, a limit is set at 
m Vlt < 170 keV. 107 - 108 

The mass for the v T is obtained from the kinematics of r decays. The r typically 
decays to many hadrons. However, the four vectors for each of the hadrons can be 
summed. Then the decay can be treated as a two-body problem with the neutrino as 
one 4- vector and the sum of the hadrons as the other vector. At this point, the same 
method described for the can be applied. Measurements are again error-limited, 
so a limit on the mass is placed. The best limit, which is m Vr < 18.2 MeV, comes 
from fits to r~ — > 27r _ 7T + ^ T and r~ — > 37r~27r + (7r°)z/ T decays observed by the ALEPH 
experiment. 109 

The experimental situation for the v e mass measurement is more complicated. The 
endpoint of the electron energy spectrum from tritium (3 decay is used to determine 
the mass. Just as in the case of the and v T , the experiments measure a value of m 2 . 
The problem is that the measurements have been systematically negative. A review 
of measurements, as a function of time, is given in Tab. 4. Recent measurements 
at Troitsk 111 and Mainz 112 are negative, but in agreement with zero. Following the 
Particle Data Group prescription for setting limit in the case of an unphysical results, 
m 2 = is assumed, with the quoted errors. Based on these results, one can extract a 
limit of approximately < 2 eV for the mass of the v e . 

The next big step in the measurement of neutrino mass from decay kinematics will 
come from the Katrin Experiment. 118 Katrin will use tritium beta decay to measure 
the mass of the neutrino to 0.2 eV. This does not reach the range of 0.05 eV, which 
our simplistic argument presented at the top of this section indicated. However, that 
argument assumed that the lightest neutrino had zero mass. A small offset from zero 
easily boosts the spectrum into the range observable by Katrin. On the other hand, 
Katrin is sensitive to only electron flavor. Thus, its sensitivity depends up on the 
amount of mixing of u e within the heaviest neutrino. 

Another method for measuring neutrino mass from simple kinematics is to use time 
of flight for neutrinos from supernovae. Neutrinos carry away ~ 99% of the energy from 
a supernova. The mass limit is obtained from the spread in the propagation times of 
the neutrinos. The propagation time for a single neutrino is given by 

o 

(Tfl \ 
1 + ^2 ) (^) 

where to is the time required for light to reach Earth from the supernova. Because the 
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neutrinos escape from a supernova before the photons, we do not know t em it- But we 
can obtain the time difference between 2 events: 



using the assumption that all neutrinos are emitted at the same time, one can obtain 
a mass limit of ~ 30 eV from the ~ 20 events observed from SN1987a at 2 sites. 119 ' 120 

This is actually an oversimplified argument. The models for neutrino emission are 
actually quite complicated. The pulse of neutrinos has a prompt peak followed by a 
broader secondary peak with a long tail distributed over an interval which can be 4 s 
or more. The prompt peak is from "neutronization" and is mainly v ei while all three 
neutrino flavors populate the secondary peak. However, the rate of v e escape is slower 
compared to and v T produced at the same time, because the z^ e s can experience CC 
interactions, while the kinematic suppression from the charged lepton mass prevents 
this for the other flavors. However, when all of the aspects of the modeling are put 
together, the bottom line remains the same: it will be possible to set stringent mass 
limits if we observe neutrinos from nearby supernovae. 

Some argue that cosmology provides a "direct measurement." Cosmological fits 
have sensitivity to neutrino masses, but the results are dependent on the cosmological 
parameters 121 and the model for relic neutrino production. There are many exam- 
ples of models with low relic neutrino densities which would significantly change the 
present interpretation of the cosmological data. 122 In the opinion of the author, until 
these issues are settled, cosmological measurements cannot convincingly compete with 
kinematic decays and supernova measurements, despite aggressive claims. 



Now that we know that neutrinos have mass, and thus are outside of expectations, the 
obvious question is: "what other Beyond Standard Model properties do they possess?" 
The APS Study on the Future of Neutrino Physics focussed on this question. The plan 
for attack was divided into three fronts: (1) Neutrinos and the New Paradigm, (2) 
Neutrinos and the Unexpected and 3) Neutrinos and the Cosmos. The remainder of 
this paper follows this structure. 

The consequences of the discovery of neutrino mass leads to a rich array of ideas. 
It is was beyond the scope of these lectures to cover the entire spectrum. So, in each of 
the three areas, two topics are chosen for extensive discussion. The reader is referred 
to the study 1 and the accompanying theory white paper 2 for further ideas. 

3.1 Neutrinos and the New Paradigm 

The first step in creating a "New Standard Model" is to incorporate neutrino mass. 
The simplest method is to introduce a Dirac mass, by analogy with the electron. This 
allows us to introduce a small neutrino mass, simply by arguing that the coupling to 
the Higgs is remarkably small. However, the unlikely smallness of the coupling has 
pushed theorists to look for other approaches. Among the oldest of these ideas is that 
neutrinos may be "Majorana particles," i.e., they are their own antiparticle. This leads 
to a new type of mass term in the Lagrangian. Through the "see-saw" mechanism, 




(23) 



3 Neutrinos We Would Like to Meet 
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which fits well with Grand Unified Theories, this can also give a motivation for the 
apparently small value of the neutrino masses. 

A direct consequence of the Majorana See-Saw Model is a heavy neutrino, with mass 
near the GUT scale. Because the heavy neutrino gets its mass through the Majorana 
rather than Dirac term of the Lagrangian, this neutrino was massive during the earliest 
periods of the universe, before the electroweak phase transition. The decays of such a 
heavy lepton could be CP violating. This would provide a mechanism for producing 
the observed matter-antimatter imbalance seen today. 

The tidiness of the the above theoretical ideas has caused this paradigm to emerge 
as the consensus favorite for the "New Standard Model." However, there is absolutely 
no experimental evidence for this theory at this time. We have no evidence for the 
Majorana nature of neutrinos. Nor do we have any evidence for CP violation in the 
neutrino system. The great challenge of the next few years, then, is to find any sign at 
all that this theory is correct. 

This section reviews how one introduces mass into the Lagrangian. The search for 
evidence of the Majorana nature of neutrinos though neutrinoless double beta decay 
is considered. Then, the prospects for finding evidence for CP violation is considered. 

3.1.1 How Neutrinos Might Get Their Mass 

The simplest assumption is that the neutrino mass should appear in the Lagrangian 
in the same way as for the charged fermions - via a Dirac mass term. In general, the 
Dirac mass term in the Lagrangian will be of the form 

m(4>LipR + ipRtpL)- (24) 

From the arguments presented in eqs. 1 through 6, we saw that the scalar "mass" term 
mixes the RH and LH states of the fermion. If the fermion has only one chirality, then 
the Dirac mass term will automatically vanish. For this reason, a standard Dirac mass 
term for the neutrino will require the RH neutrino and LH antineutrino states. 

To motivate the mass term, the most straightforward approach is to use the Higgs 
mechanism, as was done for the electron in the Standard Model. In the case of the 
electron, when we introduce a spin-0 Higgs doublet h + ), into the Lagrangian, we 
find terms like: 

g4e R (^u L {h + ) ] + Adh )^ + h.c, (25) 

where g e is the coupling constant and "h.c." is the Hermetian conjugate. The piece of 
this term proportional to ip eR '4 ) e L (h ) Ji , combined with its Hermetian conjugate, can be 
identified with the Dirac mass term, m e 'lp e ip e . We set (h°) = v/\^2, so that we obtain 
g(h°)tf) e if) e and m e = g e vj\[2. This is the Standard Model method for conveniently 
converting the ad hoc electron mass, m e , into an ad hoc coupling to the Higgs, g e and 
a vacuum expectation value (VEV) for the Higgs, v. Following the same procedure for 
neutrinos allows us to identify the Dirac mass term with m u = g v vj\[2. The VEV, v, 
has to be the same as for all other leptons. Therefore, the small mass must come from 
a very small coupling, g v . This implies that g e >5x 10 4 ^. 

There are several troublesome features to this procedure. The first issue which is 
often raised is: 

• Why would the Higgs coupling vary across eleven orders of magnitude (the ap- 
proximate ratio of the neutrino mass to the top quark mass)? 
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In fact, this question is rather odd. Disregarding the neutrinos, the masses of the 
charged fermions varies across six orders of magnitude (from the electron mass to 
the top mass). If six orders of magnitude do not bother anyone, why should eleven? 
Turning this around, if the Higgs couplings alreay seemed stretched in the charged 
fermion case, the neutrinos stretch the argument much further. This leads to the 
second troublesome issue, 

• Physically, what is occurring? 

The Higgs mechanism really gives little physical insight. While it does introduce mass, 
it has simply shifted the arbitrariness of the magnitude of the mass into an arbitrary 
coupling to a new field. 

These two questions have led theorists to look at other explanations for small neu- 
trino mass. It has been noted that neutrinos have the unique feature of carrying no 
electric or strong charge. Thus, neutrinos, alone among the Standard Model fermions, 
may be their own antiparticle, i.e. they may be Majorana particles. The nice conse- 
quence of this is a somewhat more motivated theory of mass for neutrinos. 

To understand this, first consider what is meant to be a Dirac versus a Majorana 
particle. If neutrinos are Dirac particles, then the v and the v are distinct particles, 
just as the electron and positron are distinct. The particle, v has lepton number +1 
and the antiparticle, v has lepton number —1. Lepton number is conserved in an 
interaction. Thus, using the muon family as an example, w& [L = +1) must produce 
fi~ (L = +1) and vs (L = — 1) must produce fi + (L = — 1). The alternative viewpoint 
is that the v and u are two helicity states of the same "Majorana" particle, which 
we can call " l / naj '." The ir + decay produces the left-handed v ma i and the ir~ decay 
produces the right-handed v ma3 . This model explains all of the data without invoking 
lepton number and has the nice feature of economy of total particles and quantum 
numbers, but it renders the neutrino different from all other Standard Model fermions. 

Saying that the neutrino is its own antiparticle is equivalent to saying that the 
neutrino is its own charge conjugate, ip c = ip. The operators which appear in the La- 
grangian for the neutrino in this case are the set (V'l, ipR, ip L , ip R ) and (V>l, V\r, V' c L) V ;C _r)- 
Certain bilinear combinations of these in the Lagrangian can be identified as Dirac 
masses (i.e. m^LipR + •••))• However, we also get a set of terms of the form: 



These are the "Majorana mass terms," which mix the pair of charge-conjugate states 
of the fermion. If the particle is not its own charge conjugate, then these terms au- 
tomatically vanish and we are left with only the Dirac terms. Dirac particles have no 
Majorana mass terms, but Majorana particles will have Dirac mass terms. 
The mass terms of the Lagrangian can be written in matrix form: 



The Dirac mass, m, is on the off-diagonal elements, while the Majorana mass constants, 
Ml, Mr are on the diagonal. To obtain the physical masses, one diagonalizes the 
matrix. 

One can now invoke "see-saw models" which motivate small observable neutrino 
masses. It turns out that GUT's motivate mass matrices that look like: 124 



(M l /2)(^ l c i/j l ) + (M R /2)(if R c ^ R ) + ... 



(26) 




(27) 




(28) 
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with m u <C M. When you diagonalize this matrix to obtain the physical masses, this 
results in two states which can be measured experimentally: 



m H ght ~ rnl/M, 

m-heavy ~ M 



(29) 
(30) 



Grand Unified Theories favor very large masses for the "heavy neutrino" (often called 
a "neutral heavy lepton"). It is argued that it is most "natural" to have M be at the 
GUT scale. If M ~ 10 25 eU, and mu g ht < 1 eV, as observed, then m u ~ 10 12 eV, or 
is at the TeV scale. This is rather high compared to masses of other leptons, but not 
so far beyond the top quark mass to regard the connection as crazy. So while some 
arbitrariness remains in this model, nevertheless there is a general feeling in the theory 
community that this is an improvement. 

In this theory neutrinos have only approximate handedness, where the light neutrino 
is mostly LH with a very small admixture of RH and the neutral heavy lepton is 
essentially RH. Thus we have a LH neutrino which is light, which matches observations, 
and a RH neutrino which is not yet observed because it is far too massive. 

3.1.2 Majorana vs. Dirac? 

How can we experimentally tell the difference between the Dirac (is, v) and Majorana 
(i> ma:) ) scenarios? One can imagine a straight-forward thought experiment. First, 
produce left-handed neutrinos in ir + decays. These may be ^s or they may be v™hS- 
Next, run the neutrino through a magic helicity-flipping device. If the neutrinos are 
Majorana, then what comes out of the flipping-device will be ■ These particles will 
behave like antineutrinos when they interact, showing the expected RH ^-dependence 
for the cross section. But if the initial neutrino beam is Dirac, then what comes out of 
the flipping-device will be right-handed vs, which are sterile. They do not interact at 
all. Such a helicity-flipping experiment is presently essentially impossible to implement. 
If neutrinos do have mass, then they may have an extremely tiny magnetic moment 
and a very intense magnetic field could flip their helicity. But the design requirements 
of such an experiment are far beyond our capability at the moment. Therefore, at the 
moment, we do not know if neutrinos are Majorana or Dirac in nature. 

Instead, experimentalists are pursing a different route. The Majorana nature of the 
neutrino can lead to an effect called neutrinoless double (3 decay: (Z, A) — > (Z + 2, A) + 
(e~e _ ). This is a beyond-the-Standard Model analogue to double (3 decay: (Z,A) — ► 
(Z+2, A) + (e~ e~v e i> e ). Double (3 decay is a standard nuclear decay process with a very 
low rate because there is a suppression proportional to (Gf cos 9c) 4 - Therefore, in most 
cases, if the weak decay is possible, single j3 decay ((Z, A) — ► [Z + 1, A) + e~ + u e ) will 
dominate. However, there are 13 nuclei, including 136 Xe -► 136 Ba and 76 Ge -^ 76 Se, 
for which single (3 decay is energetically disallowed. In these cases double (3 decay 
with two neutrinos has been observed. 123 If the neutrino were its own antiparticle, 
then the neutrinos produced in the double (3 decay process could annihilate, yielding 
neutrinoless double (3 decay. 

If there are Majorana neutrinos, then the amplitude for 0vf3f3 is proportional to the 
square of 




(31) 
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Figure 16: Spectrum for two- neutrino double (3 decay and expected peak for neutrinoless 
double (3 decay. 

This should be contrasted with eq. 21. The 0uj3j3 searches are probing different effective 
masses than the direct searches and the two results yield complementary information. 
Like the direct searches, the possibility of seeing Qvf3f3 depends on the amount of 
electron-flavor mixed in the most massive neutrino state. If this is small, then the 
rate of decay will be very low. Thus the hierarchy of the neutrino states affects our 
ability to observe 0v(3(3. To completely untangle the Dirac vs. Majorana question, 
three different experiments - direct mass measurement, hierarchy measurement and 
Qv(3(3 measurement - may be required. 2 

Extracting mQ U gs from a measured half-life leads to a theoretical error from the 
nuclear matrix element calculations. A favored style of calculation uses the "QRPA" 
(Quasiparticle Random Phase Approximation) 125-128 model. Using 100 Mo as an exam- 
ple, different matrix elements from QRPA calculations cause iJiQ^gg to vary by up to 2 
eV for a half-life of 4.5 x 10 23 years. 129 So the error is significant. 

The 0vf3f3 events must be separated from the standard two-neutrino double (3 
{2uf3[3) decay background. This can be done through simple kinematics cuts. The 
two-body nature of 0v[3[3 decay will cause a peak at the endpoint of the 2i/j3j3 decay 
(4-body) spectrum, as shown in Fig. 16. An advantage of observing 2v(3(3, however, 
is that measurement of its half-life allows direct measurement of the matrix element. 
At this point the 1vf3f3 decay spectrum has been observed in 10 elements. In some 
cases, such as 100 Mo, the the 2v(3(3 half-life is well measured and can be used to con- 
strain nuclear matrix element calculation. For this case, NEMO-3 reports a half life of 
(7.68 ± 0.02(stat) ± 0.54(sys)) x 10 18 y. 130 

At present, no signal for Ou/3/3 decay has been clearly observed. The present 90% 
CL limit on the lifetime from CUORICINO on 130 Te is 1.8 x 10 24 years, corresponding 
to limit of movpp < 0.2 - 1.1 eV. 131 The NEMO-3 experiment has set 90% CL limits of 
4.6 x 10 23 and 1.0 x 10 23 on 100 Mo and 82 Se, respectively. 129 The corresponding limits 
on mov/3/3 are 0.7-2.8 eV and 1.7-4.9 eV. 129 There is a candidate signal observed at 4.2<r 
from a Germanium detector, 132 although the statistical significance is under debate. 2 
The measured half-life was 1.19 x 10 25 years. Until this result is confirmed by further 
experiments, it is best to reserve judgment. 
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Figure 17: Example of two diagrams for Neutral Heavy Lepton decay which can interfere to 
produce CP violation. 

Luckily, a range of future Qv{3/3 decay experiments are on the horizon. These 
are expected to probe an order of magnitude further in lifetimes. In particular, the 
germanium-based GERDA experiment, 133 will turn on soon and will address the ex- 
istence of the possible signal. CUORE, 134 SuperNEMO, 135 EXO, 136 Majorana 137 and 
Moon 138 will extend the search even further using a wide range of elements. The reach 
of these near future 0uj3j3 covers the prediction for the inverted mass hierarchy. 

3.1.3 CP Violation in the Neutrino Sector 

An intriguing aspect of the "new Standard Model" is the heavy GUT-scale neutrinos 
which gain mass through the Majorana terms in the Lagrangian. There could be more 
than one, and likely, given the trend in the Standard Model, there would be three, 
so we can label these Ni, N2 and -/V3. These heavy neutrinos have mass prior to the 
electroweak phase transition in which the Dirac terms appear. As a result, prior to 
the electroweak phase transition, decays shown in Fig. 17 are possible. Both decays 
produce the same final state, N\ — ► IH, where t and H are oppositely charged. These 
diagrams interfere, and can lead to a different decay rate to l~ and £ + , which is CP 
violation. 

This form of CP violation would lead to a lepton asymmetry in the early universe 
which could be transferred into a baryon asymmetry. A mechanism for this already 
appears in the Standard Model, in which but B, baryon number, and L, lepton number 
are not conserved, but the difference, B — L, is exactly conserved. B and L violation 
occurs in transitions between vacuum states at high energies, called the "sphaleron 
process." Variations on this mechanism, called "leptogenesis," may explain the matter- 
antimatter asymmetry we see today. 
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N±, N2, and N3 are far too massive to be produced at accelerators in the near 
future. Thus observing CP violation in their decays is out of the question. However, 
observing CP violation in the light neutrino sector would be a plausible hint that the 
theory is correct. 

To incorporate CP violation into the three-light-neutrino model, the leptonic mixing 
matrix is expanded and written as: JJ mth CP = VK. In this case, V is very similar to 
the U of Eq. 12, but with a CP violating phase, 5: 

(C12C13 S12C13 si3e~* 5 \ 

-S12C23 - Cl 2 S 2 3Sl3e l<5 Cl 2 C 2 3 - Sl 2 S23Sl3e tS S23C13 J • (32) 
Sl2«23 - C 12 C23Sl3e tS -C12S23 - Si 2 C 2 3Sl3e tS C23C13 / 

This is analogous to the CKM matrix of the quark sector. The other term, 

K = diag(l,e^ 1 ,e^ 2+5) ) (33) 

has two further CP violating phases, 4>i and <p2- 

Now, we potentially have three non-zero CP violating parameters in the light neu- 
trino sector, 5, 4>i and 4>2i as well as one or more CP violating parameters in the heavy 
neutrino sector, where the number depends upon the total number of N. In the La- 
grangian, these all come from a matrix of Yukawa coupling constants. In principle, all 
of these phases can take on the full range of values, including exactly zero. However, 
it is difficult to motivate a theory in which some are nonzero and some are exactly 
zero. It is expected that these parameters will either all have non-zero values or all be 
precisely zero. If the latter case, then the difference between the lepton sector, with 
no CP violation, and quark sector, with clear CP violation, must be motivated. As 
a result, observation of CP violation in the light neutrino sector is regarded as the 
"smoking gun" to CP violation in the heavy sector. 

Returning to the light neutrino sector, how can the CP phases be measured? The 4> 
phases arise as a direct consequence of the Majorana nature of neutrinos. Therefore, in 
principle, the the (j) phase associated with the electron family is accessible in neutrino- 
less double beta decay. In practice, this will be extremely hard to measure because this 
term manifests itself as a change in the sum in eq. 31, which is proportional to the 0^/3/3 
decay amplitude. Thus one seeks to measure a deviation of the (as-yet-unmeasured) 
0v(3(3 lifetime from the prediction which depends upon the mixing angles (with rela- 
tively large errors at present), the (unknown) neutrino masses, and the (poorly known) 
nuclear matrix element. Even if the effect is large, observation of the effect is clearly 
hopeless in the near future. On the other hand, 5, the "Dirac" CP violating term in 
V may be accessible though oscillation searches. 

CP violation searches involve observing a difference in oscillation probability for 
neutrinos and antineutrinos. Only appearance experiments can observe CP viola- 
tion. A difference between oscillations of neutrinos and antineutrinos in disappearance 
searches is CPT violating. In oscillation appearance searches, the K matrix does not 
affect the oscillation probability because this diagonal matrix is multiplied by its com- 
plex conjugate. On the other hand, non-zero 5 can be observed. To test for non-zero 
5, the oscillation probability must depend upon the U e 3 component of Eq. 32. In other 
words, the search needs to involve transitions from or to electron flavor and involve 
the mass state V3. This combination of requirements - appearance signal, electron 
flavor involvement, and V3 mass state involvement - leads to one experimental option 
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at present: comparison of — ► v e to z-^j — > z/ e at the atmospheric Am 2 , which is Am 2 3 . 
The oscillation probability is given by: 

Piong-baseiine ^ sin 2 26>i 3 sin 2 6» 2 3 sin 2 A 

=F a sin 2#i3 sin <5cp cos #13 sin 26*12 sin 26*23 sin 3 A 

+ a sin 2#i3 cos 5c p cos #13 sin 2#i2 sin 2623 cos A sin 2 A 

+ a 2 cos 2 6*23 sin 2 29 12 sin 2 A, (34) 

where a and A are defined in eq. 20. The second term is negative for neutrino scattering 
and positive for antineutrino scattering. 

Unfortunately, eq. 34 convolutes two unknown parameters, the sign of Am? 3 (the 
mass hierarchy) and the value of 6*13, with the parameter, 5, that we want to measure. 
The problem of the mass hierarchy can be mitigated by the experimental design. The 
sign of Am 2 3 = m 2 — m 2 affects the terms where Am 2 is not squared. These terms 
arise from matter effects and can so be reduced if the pathlength in matter is relatively 
short. For long baseline experiments, which must shoot the beam through the Earth, 
this means that L must be relatively short. In order to retain the same L/E and, 
hence, the same sensitivity to Am 2 3 , E must be comparably reduced. On the other 
hand, the problem of #13 cannot be mitigated. From eq. 32, one sees that we are in 
the unfortunate situation of having the CP violating term multiplied by sin #13. The 
smaller this factor, the harder it will be to extract 5. If sin 2 26>i 3 is smaller than ~ 0.01 
at 90% CL, then substantial improvements in beams and detectors will be required. 

Eq. 34 also depends on two other as-yet-poorly understood parameters, 6*23 and 
the magnitude of Am| 3 . Disappearance experiments measure sin 2 26*23 = I-OO^q'}^, 67 
thus there is an ambiguity as to whether #23, is larger or smaller than 45°. Am^ is 
only known to about 10%. 56 ' 67 This measurement is extracted from the location of 
the "dip" in the rate versus L/E distribution of disappearance experiments, and is 
already systematics-dominated. Improvement requires experiments with better energy 
resolution 56 and better understanding of the CCQE and background cross sections. 139 
These errors lend a significant error to the analysis. 

Lastly, it is difficult to measure a v e signal which is at the ~ 1% level. Most 
beams have a substantial v e intrinsic contamination from /j, and K decays. Given that 
013 is small, this contamination is a serious issue. One solution to this problem is to 
go to an off-axis beam design. This relies on the tight correlation between energy and 
off-axis angle, 9, in two-body decays. For pion decay, which dominates most beams, 

0.43^ 

K = TTW 2 ' (35) 

where 7 = E n /m n is the Lorentz boost factor. The solution for two-body K decay 
replaces 0.43 with 0.96. Thus at 6 = 0, the relationship between E v and E w is linear. 
However, for larger 0, above a moderate energy threshold, all values of E w map to 
the same E v . This is illustrated in Fig. 18. The result is that an off axis beam 
which comes largely from pion decay is tightly peaked in energy while the v e intrinsic 
background is spread across a range of energies. This also helps to reduce events 
which are mis-reconstructed as v e scatters, such as NC 7r° production where a photon is 
lost. These "mis-ids" tend to be spread across a range of energies, since the true energy 
is misreconstructed. Thus a peaked signal, as one expects from an off-axis beam, is 
helpful in separating signal and background. 
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Figure 18: Neutrino energy versus angle off-axis for various values of pion energy. In this 
example one can see that for moderate off-axis angles, between 15-30 mrad, all pion energies 
between 3 and 20 GeV map to approximately 1 GeV neutrino energy. 

The major problem with this design is that off axis beams have substantially lower 
flux. The flux scales 140 as: 

F = (36) 

In this equation, A is the area of the detector and z is the distance to the detector. Two 
future long baseline experiments, NOvA 40 and T2K 141 are proposing off-axis beams 
for a f M — > v e search. Because of the low flux, very large detectors are required. 

In summary, the path to a test for non-zero 5 is clear but will take several steps and 
requires some luck. First, a clean measure #13 from Double Chooz and Daya Bay is 
needed. If sin 2 2#i3 < 0.005 at 90% CL, then a significant measurement of CP violation 
is unlikely to be possible in the near future. At the same time, improvements in the 
023 and Am 2 3 from disappearance (uu — > Vu) measurements at MINOS and T2K will 
improve the situation. T2K and NOvA 40 may be able to make a first exploration 
of CP parameter space, from v e appearance measurements, depending on statistics. 
NOvA may also be able to address the mass hierarchy question. This will open up 
the possibility of measuring CP violation to the next generation of very long baseline 
experiments. 142 The most sensitive of these use a beam originating at Fermilab and a 
LAr detector located at Ash River, Minnesota 45 or a Cerenkov or LAr detector located 
at the Deep Underground Science Laboratory at Homestake. 

At some point in the future, a beta beam or a neutrino factory beam could provide 
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an intense source of v e and v e fluxes, allowing comparison of v e — ► f M to u e — > z^. 
In this case one would search for events with wrong-sign muons in a calorimeter-style 
detector. This would be a striking signature with low background, especially in the 
case of a beta beam. This could allow a very precise measurement of J. 144 



3.2 Neutrinos and the Unexpected 

While it is nice to have a tidy, well-motivated theory of neutrino masses, it is discon- 
certing to have essentially no experimental evidence for this theory. Moreover, neutrino 
theories have a history of being incorrect. Only a decade ago, most theorists would 
have told you that neutrinos have no mass. Those who thought neutrinos might have 
mass believed it would be relatively large (> 5 eV), explaining dark matter. Most the- 
orists also believed that if the solar neutrinos were experiencing oscillations, the correct 
solution would be the small mixing angle MSW solution, because the mixing matrix 
should look like the quark matrix. Using the same logic, the atmospheric neutrino 
deficit, which could only correspond to large mixing angle, was routinely dismissed as 
an experimental effect. 

On the basis of this, it is wise not to constrain ourselves to the "New Paradigm." 
The reason the APS neutrino study chose to devote a chapter to "the Unexpected" 
was to emphasize the importance of being open to what nature is telling us about 
neutrinos. There are two ways to approach this idea: (1) the theory-driven approach: 
explore for properties which could, theoretically, exist and (2) the experiment-driven 
approach: follow up on anomalous results which have been observed. 

For lack of time, I will only briefly consider two examples of the first case: searching 
for a neutrino magnetic moment and searching for CPT violation. In the Standard 
Model, the neutrino magnetic moment is expected to be ~ 10 _19 ^£. Laboratory 
experiments and astrophysical limits are many orders of magnitude away from this 
level. 145 Nevertheless, if a new experiment could advance this measurement by an order 
of magnitude, that would be worth pursuing. A more startling discovery would be a 
difference in the oscillation disappearance probability of neutrinos versus antineutrinos. 
In a three-neutrino model, a difference in the rate of disappearance of neutrinos and 
antineutrinos would imply CPT violation. MINOS will be the next experiment to 
pursue such a search. 36 If CPT violation were discovered we would need to rethink the 
very basis of our theory. However, there are theorists exploring these ideas. 

The remainder of this section will focus on the second approach, exploring "anoma- 
lies" which have appeared in various experiments. Physicists today are always cautious 
about pursuing deviations from the Standard Model. Most do not, in the end, point 
to new physics. The Standard Model has been very resilient Most anomalies are ar- 
guably more likely due to systematic effects or statistical fluctuations, than to new 
physics. However, those which do "pan out" completely change the way we think. The 
solar neutrino deficit is a perfect example. So, if a new, unexpected result withstands 
questions by the community on the systematics of the experiment, then the anomaly 
becomes worth pursuing further. 

There are several examples of > 3a unexpected results in the neutrino sector which 
are worth pursuing and two cases are covered here. The first, the LSND anomaly, is be- 
ing actively pursued. The second, the NuTeV anomaly, will require a new experiment. 
Unlike most of the topics in these lectures, the NuTeV anomaly is not directly related 
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Figure 19: Left: LSND allowed range compared to short baseline experiment limits. Right 
Allowed range from the Karmen-LSND joint analysis. 149 



to neutrino oscillations and neutrino mass, and so expands the discussion, which has 
so far been rather narrowly focussed. 

Along with the known discrepancies which have reached the level of full-fledged 
anomalies, there are also examples of "unexpected results to watch." These results 
which have not yet reached the 3a level, but are showing interesting trends. For exam- 
ple, unconstrained fits to atmospheric oscillation data from a wide range of experiments 
consistently result in sin 2 2#23 best fit value greater than unity. While in each case, the 
best fit is ~ la from unity, it is the trend which is interesting, since the experiments 
involved are all very different. There is simply not enough space to cover this and other 
examples of "results to watch." 

The take-away message of this section is: the neutrino sector is a rich place for new 
physics to appear, and physicists need to be alert and open-minded to what nature is 
saying. 

3.2.1 The LSND Anomaly 

The LSND experiment ran at the LAMPF accelerator at Los Alamos National Lab- 
oratory between 1993 and 1998. The decay-at-rest (DAR) beam was produced by 
impinging 800 MeV protons on a beam dump. These produced ir + s which stop and 
decay to produce [i + s, which also stop and decay to produce and v e . These were 
studied in the range of 20 to 55 MeV. The tt~s capture, so the beam has a < 8 x 10 -4 
contamination of v e . The neutrino events were observed in a detector located 30 m 
downstream of the beam dump. 1220 phototubes surrounded the periphery of a cylin- 
drical detector filled with 167 tons of mineral oil, lightly doped with scintillator. The 
signature of a l> e appearance was v e + p —* e + n. This resulted in a two-component 
signature: the initial Cerenkov and scintillation light associated with the e + , followed 
later by the scintillation light from the n capture on hydrogen, producing a 2.2 MeV 
7. The experiment observed 87.9 ± 22.4 ± 6.0 events, a 4a excess 146 above expectation. 
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LSND is a short baseline experiment, with an L/E ~ 1 m/MeV. Thus, from the 
two-generation oscillation formula, eq. 10, one can see that this experiment is sensitive 
to Am? > 0.1 eV 2 . Other experiments have searched for oscillations at high Am 2 , 
and the two most relevant to LSND are Karmen 147 and Bugey. 148 The KARMEN 
experiment, which also used a DAR muon beam and was located 17.7 m from the 
beam dump, had sensitivity to address only a portion of the LSND region, and did 
not see a signal there. Since the design of the experiments are very similar, one can 
think of the Karmen experiment as a "near detector," which measures the flux before 
oscillation. The results were combined in a joint analysis performed by collaborators 
from both experiments; the allowed range for oscillations is shown in Figure 19. 149 
The Bugey experiment was a reactor-based v e disappearance search which set a limit 
on oscillations. Because this is disappearance and not explicitly v e — > v^, its limit is 
applicable to LSND in many, but not all, oscillation models. This limit is shown in 
Figure 19. 

Why can't we fit LSND into the three-neutrino theory? The LSND 
signal cannot be accommodated within the standard three-neutrino picture, given 
the solar and atmospheric oscillations. To see the incompatibility, first consider the 
case where the oscillation signals de-couple into, effectively, two-generation oscillations 
(eq. 10). For three generations, then Am^ = Am^ + Am^, which is clearly not the 
the case for these three signals. The more general case allows the atmospheric result to 
be due to a mixture of high (LSND-range) and low (solar-range) Am 2 values. In order 
for this model to succeed, the atmospheric Am 2 from a shape analysis must shift up 
from its present value of ~ 2 x 10~ 3 eV 2 and the chlorine experiment must have over- 
estimated the deficit of 7 Be solar neutrinos. However, the largest clash between data 
and this model arises from the Super-K v e events. This model requires that Super-K 
has missed a v e appearance signal of approximately the same size and shape as the 
deficit before detector smearing and cuts. Neutrino measurements are experimentally 
difficult and parameters do sometimes shift with time as systematics are better under- 
stood, but it seems unlikely that all of the above results could change sufficiently to 
accommodate LSND. 

Sterile Neutrinos as a Solution Additional neutrinos which do not interact 
via exchange of W or Z are called "sterile;" they may mix with active neutrinos, and 
thereby can be produced in neutrino oscillations. Experimental evidence of this would 
be the disappearance of the active flavor from the beam. In contrast to the GUT-scale 
sterile neutrinos we have already discussed, the sterile neutrinos which could explain 
LSND must be light (in the eV range), and this narrows the class of acceptable theories. 
Nevertheless, a number of possible explanations remain. 150 

Sterile neutrinos solve the LSND problem by adding extra mass splittings. The 
additional mass states must be mostly sterile, with only a small admixture of the active 
flavors in order to accommodate the limits on sterile neutrinos from the atmospheric 
and solar experiments. In principle, one might expect three sterile neutrinos. In 
practice, the data cannot constrain information on more than two sterile neutrinos. 
Therefore these are called "3+2" models. The method for fitting the data is described 
in reference. 151 One is fitting for two additional mass splittings, Am 2 4 and Am 2 5 . In 
the fit, the three mostly active neutrinos are approximated as degenerate. The mixing 
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Channel 


Experiment 


Lowest 
Am 2 


sin 2 29 
at high Am 2 


Best reach in 
sin 2 29 


^ M -> f e 


LSND 
KARMEN 
NOMAD 


0.03 eV 2 
0.06 eV 2 
0.4 eV 2 


> 2.5 x 10~ 3 

< 1.7 x 10~ 3 

< 1.4 x 10~ 3 


> 1.2 x 10~ 3 

< 1.0 x 10~ 3 

< 1.0 x 10~ 3 


z/ e disappearance 


Bugey 
Chooz 


0.01 eV 2 
0.0001 eV 2 


< 1.4 x 10" 1 

< 1.0 x lO" 1 


< 1.3 x 10~ 2 
< 5 x 10~ 2 


disappearance 


CCFR84 
CDHS 


6 eV 2 
0.3 eV 2 


NA 
NA 


< 2 x 10" 1 
< 5.3 x 10- 1 



Table 5: Results used in 3+2 fit. sin 2 29 limit is 90% CL 



matrix is also expanded by two rows and two columns. 

The data that drive the fits are the "short baseline" experiments that provide 
information on high Am 2 oscillations, summarized in Tab. 5. The combination of 

-► v e (LSND, 146 Karmen II 147 ), -»• u e (NOMAD 152 ), ^ disappearance (CDHS, 153 
CCFR84 154 ), and v e disappearance (Bugey, 148 CHOOZ 25 ) must all be accommodated 
within the model. A constraint for Super K z/„ disappearance is also included. None of 
the short baseline experiments except for LSND provide evidence for oscillations beyond 
3(j. However, it should be noted that CDHS has a la (statistical and systematic, 
combined) effect consistent with a high Am 2 sterile neutrino when the data are fit for 
a shape dependence, and Bugey has a la pull at Am 2 ~ leV 2 . As a result, these two 
experiments define the best fit combination of high and low Am 2 for the 3+2 model. 
However, there are acceptable solutions with a combination of low Am 2 values. The 
best fit, 151 has Am 2 4 = 0.92 eV 2 , Amf 5 =22 eV 2 , although there are combinations 
which work with two relatively low Am 2 values. A wide range of mixing angles can 
be accommodated, and the best fit has C/ e 4 = 0.121, = 0.204, U e 5 = 0.036 and 
[7 M 4 = 0.224. The other mixing angles involving the sterile states are not probed by 
the v„ disappearance, v e disappearance and — ► v e appearance experiments listed 
above. There is 30% compatibility for all other experiments and LSND. 

Introducing extra neutrinos, including sterile ones, would have cosmological impli- 
cations, compounded if the extra neutrinos have significant mass (>1 eV). However, 
there are several ways around the problem. The first is to note that while the best 
fit requires a high mass sterile neutrino, there are low-mass fits which work within the 
3+2 model. The second is to observe that there are a variety of classes of theories 
where the neutrinos do not thermalize in the early universe. 122 In this case, there is 
no conflict with the cosmological data, since the cosmological neutrino abundance is 
substantially reduced. 

If more than one Am 2 contributes to an oscillation appearance signal, then the data 
can be sensitive to a CP-violating phase in the mixing matrix. Experimentally, for this 
to occur, the Am 2 values must be within less than about two orders of magnitude of 
one another. 
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In 3+2 CP-violating models: 155 



pfi^W) = A\U ei \ 2 \U^\ 2 sm 2 x Al (37) 
+4|[/ e5 | 2 |?7 M 5| 2 sin 2 x 5 i 

+8|£4 4 |IM|£ / e5||^5| 

sin X41 sin X51 cos(x54 + ^54) > 



where in the last line, the negative sign is for neutrino oscillations and the positive sign 
is for antineutrino oscillations, and defined: 

= 1.27Am 2 i L/£', 54 = argiU^U^U^U*^. 

Thus the oscillation probability is affected by CP violation through the term ^54. The 
CP conserving cases are ^54 = and 180 degrees. 

MiniBooNE First Results The main purpose of the MiniBooNE experiment 
was to resolve the question of the LSND signal. First results of this experiment, 
presented in April, 2007, considered those explanations with a high expectation for 
Vji, — ► v e oscillations. This includes the CP conserving 3+2 model and many cases of 
CP violating 3+2 models described above. As will be described below, the first results 
are incompatible with — ► u e oscillations, but show an unexpected low energy excess, 
very much in keeping with the subsection title of "Neutrinos and the Unexpected." 

The MiniBooNE experiment uses the Fermilab Booster Neutrino Beam, which is 
produced from 8 GeV protons incident on a beryllium target located within a mag- 
netic focusing horn. The current of the horn can be reversed such that the beam is 
dominantly neutrinos or antineutrinos. The first results are from neutrino running. 
The MiniBooNE detector is located L = 541 m from the primary target, and the neu- 
trino flux has average energy of ~ 0.75 GeV. The detector is located 541 m from the 
front of the beryllium target and consists of a spherical tank of radius 610 cm that is 
covered on the inside by 1520 8-inch photomultiplier tubes and filled with 800 tons of 
pure mineral oil (CH2). Neutrino events in the detector produce both Cerenkov and 
scintillation light. 

In order to test the LSND result, the MiniBooNE design maintains L/E ~ 1 m/MeV 
while substantially changing the systematic errors associated with the experiment. This 
is accomplished by increasing both L and E by an order of magnitude from the LSND 
design. This changes the source of the neutrinos [y^ from energetic pions rather than 
i?n from stopped muons), the signature for the signal, and the major backgrounds in 
the detector. In its first run, in neutrino mode, MiniBooNE collected over a million 
clean, neutrino events. About 99.5% of the MiniBooNE neutrino events are estimated 
to be v ^-induced, while 0.5% are estimated to be due to "intrinsic" v e background in 
the beam. 

The initial MiniBooNE results were analyzed within an appearance-only, two neu- 
trino oscillation context. While the LSND signal must be a result of a more complex 
oscillation model, in most cases a — ► ^ e -like oscillation signal is predicted. After 
the complete v e event selection is applied, the total background was estimated to be 
358 ± 35 events, while 163 ± 21 signal events were expected for the LSND central 
expectation of 0.26% — ► v e transmutation. 
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The top plot of Fig. 20 shows candidate v e events as a function of reconstructed 
neutrino energy (E® E ). The vertical dashed line indicates the minimum E„ E used in 
the two-neutrino oscillation analysis. There is no significant excess of events (22+19+35 
events) for 475 < E® E < 1250 MeV; however, an excess of events (96 ± 17 ± 20 events) 
is observed below 475 MeV. In the top plot, the points show the statistical error, while 
the histogram is the expected background with systematic errors from all sources. The 
background subtracted excess as a function of E® E is shown in the bottom plot, where 
the points represent the data with total errors. Oscillation scenarios are indicated by 
the histograms. 

The low-energy excess cannot be explained by a two-neutrino oscillation model, and 
its source is under investigation. The low energy events isolated by the cuts, including 
the excess events, are single-ring and electromagnetic-like, with no unusual detection 
issues. The low energy excess events are neither consistent with the spatial nor energy 
distributions of photons coming from interactions outside of the tank. Nor are they 
consistent with the energy distribution from single photons from radiative A decays 
(A — ► N + 7). Mis-identification of tt° events is well constrained in MiniBooNE by the 
rate of reconstructed 7r° events, studied as a function of 7r° momentum. 156 This rate 
would need to be mis-measured by well over a factor of three to explain the excess, far 
outside of the error of the analysis. 

With that said, as shown in Fig. 20, the excess is not in agreement with a simple 
two-neutrino — > v e oscillation signal. This figure shows the predicted spectrum 
when the best-fit two-neutrino oscillation signal is added to the predicted background. 
The bottom panel of the figure shows background-subtracted data with the best-fit 
two-neutrino oscillation and two oscillation points from the favored LSND region. 

A single-sided raster scan to a two neutrino appearance-only oscillation model is 
used in the energy range 475 < E„ E < 3000 MeV to find the 90% CL limit correspond- 
ing to Ax 2 = xfimit ~ Xbestfit = l-^4- As shown by the top plot in Fig. 21, the LSND 
90% CL allowed region is excluded at the 90% CL. A joint analysis as a function of 
Am 2 , using a combined \ 2 of the best fit values and errors for LSND and MiniBooNE, 
excludes at 98% CL two-neutrino appearance-only oscillations as an explanation of the 
LSND anomaly. The bottom plot of Fig. 21 shows limits from the KARMEN 147 and 
Bugey 148 experiments. This is plot represents an example of the problem of apples- 
to-apples comparisons raised in sec. . The MiniBooNE and Bugey curves are 1-sided 
upper limits on sin 2 26 corresponding to Ax 2 = 1-64 - hence directly comparable - 
while the published KARMEN curve is a "Feldman Cousins" contour. 

Next Steps From the initial MiniBooNE result, one can draw two conclusions: 
(1) there is excellent agreement between data and prediction in the analysis region 
originally defined for the two- neutrino oscillation search and (2) there is a presently 
unexplained discrepancy with data lying above background at low energy. This com- 
bination of information severely limits models seeking to explain the LSND anomaly. 

Interpreting the MiniBooNE data as appearance-only and combining this result 
with other data in a 3+2 fit does not give a satisfactory result. 157 While MiniBooNE 
and LSND are compatible if CP violation is allowed in the 3+2 model, there is tension 
between these results and the disappearance experiments. This might be addressed if 
a 3+2 interpretation of the MiniBooNE result were expanded to include the possibility 
of is p disappearance and intrinsic u e disappearance. This analysis is underway by the 
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Figure 20: The top plot shows the number of candidate v e events as a function of . 
Also shown are the best-fit oscillation spectrum (dashed histogram) and the background 
contributions from and v e events. The bottom plot shows the number of events with the 
predicted background subtracted as a function of E® E . The two histograms correspond to 
LSND solutions at high and low Am 2 . 

MiniBooNE collaboration. 158 Most likely, if a good fit is obtained in a 3+2 scenario, 
it will require some level of CP violation. MiniBooNE is presently collecting data in 
antineutrino mode. However, this is a small data set (~ 2 x 10 20 protons on target 
producing the beam) and future running to reach roughly three times the statistics 
will be required to make a decisive statement. Other, alternative explanations are also 
being explored 159-161 . 

An upcoming result which will shed light on the question is the analysis of the 
MiniBooNE data from the NuMI beam. This beam is 110 mrad off-axis, with a ir peak 
of average energy of about 200 MeV and a K peak of about 2 GeV, and a length of 
750 m. If an excess of events is observed in this analysis, this rules out mis-estimate 
of intrinsic v e in the Booster Neutrino Beam as the source of the MiniBooNE excess. 
Results from this study are expected in autumn, 2007. 

If the unexplained excess persists after the above studies, then it will be valuable 
to introduce a detector which can differentiate between electrons and photons. That 
is the goal of MicroBooNE, 44 which uses a Liquid argon TPC (LArTPC) detector. 
This is particularly sensitive at low energies and nearly background-free. Specifically, 
this detector has a v e efficiency > 80% and rejects photons efficiently through dE/dx 
deposition in the first ~2 cm of the shower. With these qualities the detector can 
be an order of magnitude smaller in size than MiniBooNE, making quick construction 
feasible. A proposal for this experiment will be submitted to the Fermilab PAC in 
autumn 2007. 
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Figure 21: The top plot shows the — > v e MiniBooNE 90% CL limit (thick solid curve) 
and sensitivity (dashed curve) for events with 475 < < 3000 MeV. Also shown is the 
limit from a second cross-check analysis (thin solid curve). The bottom plot shows the limits 
from the KARMEN 147 and Bugey 148 experiments. The shaded areas show the 90% and 99% 
CL allowed regions from the LSND experiment. 
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Table 6: left and right handed coupling constants. 



3.2.2 The NuTeV Anomaly 

Neutrino scattering measurements offer a unique tool to probe the electroweak inter- 
actions of the Standard Model (SM). The NuTeV anomaly is a 3<r deviation of sin 2 9w 
from the Standard Model prediction. 162 sin 2 9w parameterizes the mixing between the 
weak interaction Z boson and the photon in electroweak theory. Deviations of measure- 
ments of this parameter, and its partner parameter, p, the relative coupling strength 
of the neutral-to-charged-current interactions, may indicate Beyond-Standard-Model 
physics. This section also highlights that fact that new neutrino properties may be 
revealed in TeV-scale interactions at LHC, which has not been addressed previously. 

The NuTeV experiment represents a departure from the previous train of thought 
in several ways. NuTeV was a deep inelastic neutrino scattering experiment, and thus 
is performed at significantly higher energy than the experiments previously discussed. 
Also, while NuTeV did an oscillation search, it was mainly designed for another pur- 
pose: precision measurement of electroweak parameters. We will focus on that purpose 
here. As a result, this analysis allows new issues related to neutrino physics to be 
brought into the discussion. 

sin 2 9w in Neutrino Scattering and Other Experiments 

In neutrino scattering, the neutral current cross section depends upon sin 2 9w- The 
dependence is a function of the neutrino flavor and the target. NuTeV was a muon- 
neutrino-flavor scattering experiment. In this case, the NC cross sections for scattering 
from a light fermion target are: 

OfM^Il _ ^M^-,)-)^*)-, (38) 

In this equation, / is the type of light fermion: / = e~ ,u,d,s,c. Several constants 
appear: Gf is the Fermi constant, Mz is the mass of the Z. The two kinematic 
variables are: s, the effective center of mass energy, which depends on the mass of / 
and y, the inelasticity (see definitions in Sec. 1.2). £f,r/ are left and right handed 
coupling constants which are given in Table 6. 

While neutrino scattering has traditionally been a method for measuring sin 2 9\y, 
the "Standard Model Prediction" quoted in literature comes from the very precise 
measurements made by the LEP and SLD experiments, which have been summarized 
by the Electroweak Working Group. 163 sin 2 9y/ appears in various measurements from 
e + e~ scattering at the Z pole. An example which leads to a highly precise measurement 
is the "left-right asymmetry" measured from polarized scattering at SLD: Am = (cjl — 
o~r) I (o-l + pr) where &l and an refer to the scattering cross sections for left- and right- 
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polarized electrons, respectively. In this case the asymmetry is given by: 



A L r(Z°) = ± r =r, 40 

(^- S in 2 ^r)) 2 +sin 4 4 ff )' 



When comparing quoted values of the weak mixing angle sm 2 6w: care must be 
taken because this parameter is defined in various ways. The simplest definition is the 
"on shell" description: 

1-M^/Ml^sin 2 ^ n - Shcll) . (41) 

This is the definition commonly used in neutrino physics. In the discussion which 
follows, if not explicitly labeled, the on-shell definition for sin 2 6y/ is used. A variation 
on this definition uses the renormalized masses at some arbitrary scale fi which is 
usually taken to be Mz- 

1 - Mwirf/Mz^) 2 = sin 2 ^~ S) . (42) 

However, the LEP experiments used the "effective" weak mixing angle which is related 
to the vector and axial vector couplings: 

J(l-^M)^sin 2 4 ff ). (43) 

This is what appears in eq. 40 above. One must convert between definitions, which 
have different radiative corrections and renormalization prescriptions, in order to make 
comparisons. 

The parameter sin 2 Oyy evolves with Q 2 , the squared 4-momentum transfer of the 
interaction. Fig. 22 illustrates this evolution. The highest Q 2 measurements are from 
LEP and SLD, with Q 2 = M§. There are several types of experiments, including 
neutrino experiments, which measure sin 2 9w with Q 2 <C m 2 z . NuTeV was performed at 
Q 2 = 1 to 140 GeV 2 , (Q 2 ) = 26 GeV 2 , (Q?) = 15 GeV 2 . The lowest Q 2 measurements 
are from studies of atomic parity violation in the nucleus 164 (APV), which arises due 
to the electroweak interference of the photon and the Z in the boson exchange between 
the electrons and the nucleus. This samples Q 2 ~ 0. At higher Q 2 , there is the result 
from SLAC E158, a M0ller scattering experiment at average Q 2 = 0.026 GeV 2 . 165 
Using the measurements at the Z-pole with Q 2 = M 2 to fix the value of sin 2 9w , an d 
evolving to low Q 2 , Fig. 22 166 shows that APV and SLAC E158 are in agreement with 
the Standard Model. 

NuTeV is strikingly off the prediction of Fig. 22. Neutrino scattering may measure 
a different result because new physics enters the neutrino process differently than the 
other experiments. Compared to the colliders, neutrino physics measures different 
combinations of couplings. Also neutrino scattering explores new physics through 
moderate space-like momentum transfer, as opposed to the time-like scattering at the 
colliders. With respect to the lower energy experiments, the radiative corrections to 
neutrino interactions allow sensitivity to high-mass particles which are complementary 
to the APV and M0ller-scattering corrections. 

The NuTeV Result The NuTeV experiment provides the most precise measure- 
ment of sin 2 6w from neutrino experiments. The measurement relied upon deep in- 
elastic scatter (DIS). It was performed using a "Paschos-Wolfenstein style" 167 analysis 
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Figure 22: Measurements of sin 2 6 W as a function of Q. 166 The curve shows the Standard 
Model expectation. 

which is designed to minimize the systematic errors which come from our understanding 
of parton distributions and masses. 

This method requires separated v and v beams. In this case, the following ratios 
could be formed: 

= ?MC (44) 
°cc 

R 9 = % (45) 



°cc 



(46) 



Paschos and Wolfenstein 167 recast these as: 



R -^h^r c -^^' <47) 

where r = ^cc/ a CC- ^ n ^he case °^ > many systematics cancel to first order. 
In particular, the quark and antiquark seas for u, d, s, and c, which are less precisely 
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known than the valence quark distributions, will cancel. Charm production only enters 
through d va ience which is Cabbibo suppressed and at high x, thus the error from the 
charm mass is greatly reduced. One can also form R + , but this will have much larger 
systematic errors, and so the strength of the NuTeV analysis lies in the measurement 
of R-. 

According to the "Paschos-Wolfenstein" method, an experiment should run in neu- 
trino and antineutrino mode, categorize the events as CC or NC DIS, and then form 
R~ to extract sin 2 #^- This requires identifying the CC or NC events properly in 
NuTeV's iron-scintillator/drift-chamber calorimeter. Most CC DIS events have an ex- 
iting muon, which causes a long string of hits in the scintillator and are therefore called 
"long." Most NC DIS events are relatively "short" hadronic showers. However, there 
are exceptions to these rules. A CC event caused by interaction of an intrinsic v e in 
the beam will appear short. An NC shower which contains a pion-decay-in-flight, pro- 
ducing a muon, may appear long. The connection between long vs. short and CC vs. 
NC must be made via Monte Carlo. 

NuTeV measurement is in agreement with past neutrino scattering results, although 
these have much larger errors. However, the NuTeV result is in disagreement with the 
global fits to the electroweak data which give a Standard Model value of sin 2 9w = 
0.2227. 162 

Explanations In the case of any anomaly, it is best to start with the commonplace 
explanations. Three explanations for the NuTeV anomaly that are "within the Stan- 
dard Model" have been proposed: the QCD-order of the analysis, isospin violation, 
and the strange sea asymmetry. The NuTeV analysis was not performed at a full NLO 
level. However, the effect of going to NLO on NuTeV can be estimated, 168 and the 
expected pull is away from the Standard Model. The NuTeV analysis assumed isospin 
symmetry, that is, u{x) v = d{x) n and d{x) p = u{x) n . Various models for isospin viola- 
tion have been studied and their pulls range from less than la away from the Standard 
Model to ~ la toward the Standard Model. 169 Variations in the strange sea can either 
pull the result toward or away from the Standard Model expectation, 169 but not by 
more than one sigma. 

With respect to Beyond-Standard-Model explanations, Chapter 14 of the APS Neu- 
trino Study White Paper on Neutrino Theory 2 is dedicated to "The Physics of NuTeV" 
and provides an excellent summary. The discussion presented here is drawn from this 
source. 

The NuTeV measurements of R u and R u , the NC-to-CC cross sections, are low 
compared to expectation. For this to be a Beyond-Standard-Model effect, it therefore 
requires introduction of new physics that suppresses the NC rate with respect to the 
CC rate. Two types of models produce this effect and remain consistent with the 
other electroweak measurements: (1) models which affect only the Z couplings, e.g., 
the introduction of a heavy Z' boson which interferes with the Standard Model Z; or 
(2) models which affect only the neutrino couplings, e.g., the introduction of moderate 
mass neutral heavy leptons which mix with the neutrino. 

Any Z' model invoked to explain NuTeV must selectively suppress NC neutrino 
scattering, without significantly affecting the other electroweak measurements. This 
rules out most models, which tend to increase the NC scattering rate. Examples of 
successful models are those where the Z' couples to B — SL^ 170 or to L M — L T . 171 
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Moderate-mass neutral heavy leptons, a.k.a. " neutrissimos," can also produce 
the desired effect. Suppression of the coupling comes from intergenerational mixing of 
heavy states, so that the is a mixture: 

v,j, = (cos a) flight + (sin a Heavy (48) 

The Zv^v^ coupling is modified by cos 2 a and the Wfiu^ coupling is modified by cos a. 
Neutrissimos may have masses as light as ~ 100 GeV. 172 These new particles can play 
the role of the seesaw right-handed neutrinos, as long as one is willing to admit large 
tuning among the neutrino Yukawa couplings. 172 So this offers an alternative to the 
GUT-mass heavy neutrino model discussed in sec. 3.1. 

If neutrissimos exist, they would be expected to show up in other precision exper- 
iments. One must avoid the constraints on mixing from Of/3/3 (recall eq. 31 to see 
why these experiments have sensitivity to the mixing) . These experiment place a limit 
of \U e 4\ 2 at less than a few xlO -5 for a 100 GeV right-handed neutrino. Rare pion 
and tau decays constrain |?7 M 4| 2 to be less than 0.004 and \U T 4\ 2 to be less than 0.006, 
respectively. 

Neutrissimos would be produced at LHC, thus neutrino physics can be done at the 
highest energy scales! However, they may be difficult to observe. One would naturally 
look for a signal of missing energy. However, neutrissiomos will not necessarily decay 
invisibly; for example one can have N — > I + W and the W may decay to either two 
jets or a neutrino-charged-lepton pair. Only the latter case has missing energy. This 
may make them difficult to identify. 

If the neutrissimo is a Majorana particle, then these could provide a clue to the 
mechanism for leptogenesis. The present models of leptogenesis require very high mass 
scales for the neutral lepton. However, theorists are identifying ways to modify the 
model to accommodate lower masses. 173 There also may be a wide mass spectrum 
for these particles, with one very heavy case that accommodates standard leptogenesis 
models, while the others have masses in the range observable at LHC. 174 

NuSOnG and Other Possibilities A new round of precision electroweak mea- 
surements can be motivated by the NuTeV anomaly as well as the imminent turn-on 
of LHC. These measurements are best done using neutrino-electron scattering, because 
this removes the quark-model related questions discussed in the previous section. Two 
possible methods for such a measurement are f„ scattering with higher statistics, using 
a NuTeV-style beam, or a u e s scattering measurement from a reactor. In either case, to 
provide a competitive measurement, the error from the best present neutrino-electron 
scattering measurement, from CHARM II, must be reduced by a factor of five. 

The NuSOnG (Neutrino Scattering On Glass) Experiment 175 is proposed to run 
using a beam produced by 800 GeV protons on target from the TeVatron. The plan 
is to use a design which is inspired by the CHARM II experiment: a target of SIO2 
in one quarter radiation length panels, with proportional tubes or scintillator to allow 
event reconstruction. The detector will have a 2.6 kton fiducial volume. The major 
technical challenge of such an experiment is in achieving the required rates from the 
TeVatron, as x20 the rate of NuTeV proton delivery is required. 

Alternatively, a measurement of the weak mixing angle using anti-neutrinos from 
reactors may be possible. 176 The weak mixing angle can be extracted from the purely 
leptonic v e e "elastic scatter" (ES) rate, which is normalized using the T> e p "inverse 
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beta decay" (IBD) events, to reduce the error on the flux. Thus, a hydrocarbon 
(scintillator oil) based detector, which has free proton targets for the IBD events, is 
ideal. Gadolinium (Gd) doping is necessary for a high rate of neutron capture, which 
constitutes the signal for the IBD events. A window in visible energy of 3 to 5 MeV is 
selected to reduce backgrounds from contaminants in the oil and cosmic-muon-induced 
isotopes. In this energy range, the dominant contamination comes from the progeny of 
the uranium and thorium chain. This would clearly be an ambitious, state-of-the-art 
measurement, but could be done at a new reactor experiment where the detector is in 
close proximity to the source and had high shielding from cosmic rays. 

3.3 Neutrinos and the Cosmos 

Neutrinos are ubiquitous in the universe, and their presence and interactions must be 
incorporated into astrophysical and cosmological models. Nearly any new neutrino 
property will have direct consequences in these fields, which must be examined. As 
an illustrative example, the first discussion considers the impact of introducing of 
relatively light sterile neutrinos to the theory. Sterile neutrinos with keV-scale masses 
can explain dark matter as well other astrophysical questions. 

As we improve our capability for detecting astrophysical neutrinos, these become a 
new source of neutrinos for study. The second example is a case in point: the search 
for ultra- high energy sources of neutrinos. The discovery of such sources would be of 
great interest to astrophysics, and the particle physics we can do with such a "beam" 
is remarkable. 

While this section concentrates on the unknown, it is interesting to note that the 
known astrophysical sources of neutrinos are sufficiently intense that these neutrinos 
are already a possible background to other physics measurements. An example is the 
case of dark matter searches, which aim to measure cross sections as small as 10~ 46 cm 2 . 
These experiments will have to contend with the background from coherent scattering 
(y + N — ► v + N) of solar neutrinos, which has a cross section of 10~ 39 cm 2 and which 
produces a recoil nucleon that is very much like the expected dark matter signal. 177 

3.3.1 Neutrinos as Dark Matter 

From the mid-1980's through mid-90's a 5 eV v T was considered a likely candidate 
for dark matter. This was the motivation for the Chorus and NOMAD search for 
Vji, — ► v T oscillations in the > 10 eV 2 range 178 ' 179 as well as the proposed COSMOS 
experiment. 180 

In the late 90's and early 2000's, two measurements led to a shift in opinion about 
neutrinos as candidates for dark matter. The first was the Super-K confirmation of 

oscillations. As discussed in sec. 2.1.4, the cleanest explanation, which fits within a 
three-neutrino model, is that this effect is — ► v T . Combining this information with 
the direct limit on the v e implied that neutrinos were unlikely to have masses in the 5 
to 10 eV range, as required for dark matter. Also, at the same time, studies of the large 
scale structure of the universe indicated that dark matter must be non-relativistic, or 
"cold." Relativistic, or "hot dark matter," like neutrinos, would smooth the large scale 
structure far beyond observations. 181 

The idea of neutrinos as dark matter candidates fell into disfavor. For some time, 
the more likely solution was assumed to be WIMPs, Weakly Interacting Massive Par- 
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Figure 23: Bounds for u e — > v s oscillations from astrophysics and cosmology. Allowed regions 
for neutrino cold, warm and hot dark matter are shown. 187 

tides. The "weak" in this name is somewhat confusing, since it does not refer strictly 
to the weak interaction - other Beyond-Standard-Model interactions are involved. It 
is simply meant to say that the interaction rate is very low. 

For some time, the lightest supersymmetric particle has been the most favored 
candidate for the WIMP. However, this is now starting to be questioned, as no evidence 
for supersymmetry has been observed at colliders. 182 This makes the formulation of the 
theory more awkward, and SUSY explanations for dark matter have been pushed from 
the "Minimal Super Symmetric Model" to the "Next-to-Minimal Super Symmetric 
Model," and even this is challenged. 183 If supersymmetry does not show up at LHC, 
then a new explanation for dark matter must be found. 

As a result, neutrino models are being reconsidered. 184-193 The new models involve 
neutrinos which are mostly sterile, with a very tiny mixing with the light flavors and 
with keV masses (0.5 keV < m v < 15keV and sin 2 26 > 10~ 12 ) Because of the high 
mass, they are not relativistic and are regarded as "warm" or "cold." At high mass, the 
large scale structure limits are much less stringent. A recent analysis from the Sloan 
Digital Sky Survey (SDSS) finds that a sterile neutrino mass above 9 keV can reproduce 
the power spectrum. 194 Only tiny mixing with the active neutrinos is required in order 
to produce the dark matter. With such small mixings, this model easily evades all 
accelerator-based bounds on v e — > v s oscillations. As shown in Fig. 23, this model also 
escapes the cosmological limits on v e — > v s from the CMB measurements, from Big 
Bang Nucleosynthesis (BBN), and supernova limits (SN), assuming negligible lepton 
number asymmetry, L, in the early universe. 

Currently the only constraints on keV sterile neutrinos come from X-ray astron- 
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omy 93 which are searching for evidence of radiative decay of the massive neutrino into 
a lighter state, v-i — > v\ + 7. This proceeds through loop diagrams where the photon is 
coupling to a If or a charged lepton in the loop. 195 Because this is a 2-body decay, one 
is searching for a spectral line in the x-ray region. For a Dirac-type sterile neutrinos of 
mass m s the decay rate is given by: 

r 7 (m s ) = 1.36 x 10- 2 V 1 (^)( I ^) 5 , (49) 

which is clearly tiny, even for a keV scale neutrino. This is important as the dark matter 
neutrinos must be stable on the scale of the lifetime of the universe. No signal has been 
observed and the current mass limit, from the Chandra X-ray telescope observations, 
ranges from >3 to >6 keV depending on model assumptions. 

Having motivated a keV-mass sterile state using dark matter, one can explore the 
consequences in other areas of cosmology and astrophysics. The small mixing allows 
these neutrinos to evade bounds from big bang nucleosynthesis. 196 Their presence may 
be beneficial to models of supernova explosions and pulsar kicks, as discussed below. 

The existence of these neutrinos may improve the supernova models substantially. 
The problem faced by most models is that the supernova stalls and fails to explode. 
As modest increase in neutrino luminosity during the epoch when the stalled bounce 
shock is being reheated (i p b < 1 s) will incite the explosion. 197 As the supernova occurs, 
neutrinos will oscillate and even be affected by MSW resonances. If neutrinos oscillate 
to a sterile state, then their transport-mean-free-paths become larger. This increases 
the neutrino luminosity at the neutrino sphere and makes "the difference between a 
dud and an explosion." 197 

Also, sterile neutrinos in the 1 to 20 keV mass range can also be used to explain 
the origin of pulsar motion. Pulsars are known to have large velocities, from 100 
to 1600 km/s. This is a much higher velocity than an ordinary star which typically 
has 30 km/s. Pulsars also have very high angular velocities. Apparently, there is some 
mechanism to give pulsars a substantial "kick" at birth, which sends them off with high 
translational and rotational velocities. One explanation for the kick is an asymmetric 
neutrino emission of sterile neutrinos during or moments after the explosion which 
forms the pulsar. 196 An asymmetry in neutrino emission occurs because the "urea 
reactions" 

u e + n <-> p + e~ , (50) 
v e +p <-> n + e + , (51) 

are affected by magnetic fields which trap electrons and positrons. If the neutrinos 
oscillate to a sterile state, they stream out of the pulsar. If the sterile neutrinos 
have high mass, they can provide a significant kick. There are a number of solutions, 
either with standard neutrino oscillations, as described by eq. 10 or with an MSW-type 
resonance. All require a sterile state in the 1 to 20 keV range with very small mixing, 
compatible with the dark matter scenario described above. 

In summary, this is an example of how introducing a new neutrino property, i.e. 
sterile companions to the known neutrinos, can have a major impact on astrophysical 
models. This is an interesting case in point, because it is unlikely that these neutrinos 
will be observable in particle physics experiments in the near future. At present, the 
only detection method is through X-ray emission due to the radiative decay. Thus this 
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is, at the moment, an example of a neutrino property which is entirely motivated and 
explored in the context of astrophysics. 



3.3.2 Ultra High Energy Neutrinos 

All of the experiments so-far discussed have used neutrino sources in the energy range 
of a few MeV to many GeV. We do not know how to produce neutrino beams at 
higher energies. However, nature clearly has high energy acceleration mechanisms, 
because cosmic rays with energies of 10 8 GeV have been measured. A new generation 
of neutrino experiments is now looking for neutrinos at these energies and beyond. 
These include AMANDA, 34 ICEcube, 198 Antares, 35 and Anita. 199 

These experiments make use of the fact that the Earth is opaque to ultra-high 
energy neutrinos. The apparent weakness of the weak interaction, which is due to the 
suppression by the mass of the W in the propagator term, is reduced as the neutrino 
energy increases. Amazingly, when you reach neutrino energies of 10 17 eV, the Earth 
becomes opaque to neutrinos. To see this, recalling the kinematic variables defined 
in sec. 1.2, consider the following back-of-the-envelope calculation. For a 10 8 GeV u, 
s = 2ME y = 2 x 10 8 GeV 2 . Most interactions occur at low x; and at these energies 
^typical ~ 0.001. For neutrino interactions, the average y is 0.5. Therefore, using 
Q 2 = sxy, we find Q 2 typical = (2 x 10 8 )(1 x 10~ 3 )(0.5) = 1 x 10 5 Gev 2 . The propagator 
term goes as 



which is approximately 10 3 for our "typical" case. The typical cross section atypical 
is: 

E x {utot/E) x (prop term) = 10 8 (GeV)(0.6 x 10- 38 cm 2 /GeV)10- 3 
= 0.6 x 10~ 33 cm 2 . From this we can extract the interaction length, Ao, on iron, by 
scaling from hadronic interactions, which tells us that at 30 mb (= 0.3 x 10~ 25 cm 2 ), 
lAo ~ 10 cm. This implies that Ao for our very high energy vs is ~ 5 x 10 3 km. 
However, the Earth is a few xlO 4 km. Thus all of the neutrinos interact; the Earth is 
opaque to them. 

This opens up the opportunity to instrument the Earth and use it as a neutrino 
target. One option is to choose a transparent region of the Earth, ice or water, and 
instrument it like a traditional neutrino detector. This has been the design chosen by 
AMANDA, ICECube and Antares, which use phototubes to sense the Cerenkov light 
produced when charged particles from neutrino interactions traverse the material. The 
largest of these detectors are on the order of (1 km) 3 of instrumented area. The 
second method exploits the Askaryan effect 200 in electromagnetic showers. Electron 
and positron scattering in matter have different cross sections. As the electromangetic 
shower develops, this difference leads to a negative charge asymmetry, inducing strong 
strong coherent Cerenkov radiation in the radio range. The pulse has unique and easy- 
to-distinguish broadband (0.2 - 1.1 GHz) spectral and polarization properties which 
can be received by detecting antennas launched above the target area. In ice, the 
radio attenuation length is 1 km. The Anita Experiment, which uses such a detector, 
can view 2 x 10 6 km 3 of volume. This makes it by far the world's largest tonnage 
experiment. 




(52) 
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Neutrinos with energies above 10 4 GeV have yet to be observed. However, they 
are expected to accompany ultra high energy cosmic rays, which have been observed. 
Nearly all potential sources of ultra-high energy cosmic rays are predicted to produce 
protons, neutrinos, and gamma rays at roughly comparable levels. Ultra high energy 
protons, which have been observed by the HiRes 201 and Auger 202 experiments are guar- 
anteed sources of neutrinos through the Gresein-Zatsepi-Kuzmin (GZK) interaction. In 
this effect, protons above Eqzk = 6 x 10 10 GeV scatter from the cosmic microwave 
background: p~f — ► A + — ► nir + . This degrades the energy of protons above Egzk, 
leading to an apparent cutoff in the flux called the "GZK cutoff." There are several 
A resonances and the CMB photons have an energy distribution, so the cutoff is not 
sharp. But it has been clearly observed by both HiRes and Auger. 203 As a result, ultra- 
high energy pions are produced and these must decay to ultra-high-energy neutrinos. 
There may be other, more exotic mechanisms for producing an ultra-high energy flux, 
possibly with energies beyond the GZK cutoff. Since, unlike the protons, these neu- 
trinos do not interact with the cosmic microwave background, they can traverse long 
distances and can be messengers of distant point sources. 

There many reviews of the exotic physics one can do with ultra high energy neutrino 
interactions. The opportunities include 204 gravitational lensing of neutrinos, the search 
for bumps or steps in the NC/CC ratio, the influence of new physics on neutrino 
oscillations at high energies, the search for neutrino decays, neutrino interaction with 
dark matter WIMPs, and the annihilation of the ultra high energy neutrinos by the 
cosmic neutrino background. This final example is interesting because significant limits 
have been set by Anita-lite, a small prototype for Anita that flew only 18.4 days. This 
illustrates the power of even a small experiment entering an unexplored frontier of 
particle physics. 

The 1 eV mass neutrino implied by the LSND anomaly could be a candidate for 
the source of the ultra-high energy cosmic rays observed on Earth. 205 This neutrino, if 
produced at ultra-high energies, could annihilate on the cosmic neutrino background 
producing a "Z-burst" of ultra high energy hadrons. This was offered as an explana- 
tion for ultra-high energy cosmic which were observed by the AGASA experiment. 206 
Scaling from the AGASA rate, a prediction for the flux of ultra-high neutrinos in the 
energy range of 10 18 5 < E v < 10 23 5 eV for Z-burst models was made. 207 ' 208 Based 
on this flux, Anita-lite was predicted to see between 5 and 50 events at >99% CL. 
During its short run, this prototype detector observed no events in the energy range 
and therefore could definitively rule out this model. 209 Shortly thereafter, the AGASA 
events were shown to be due to energy miscalibration. 210 

4 Conclusions 

The goal of this review was to sketch out the present questions in neutrino physics, and 
discuss the experiments that can address them. Along the way, I have highlighted the 
experimental techniques and challenges. I have also tried to briefly touch on techno- 
logical advances expected in the near future. This text followed the structure of the set 
of lectures entitled "Neutrino Experiments," given at the 2006 TASI Summer School. 

Neutrino physics is an amalgam of astrophysics, cosmology, nuclear physics, and 
particle physics, making the field diverse and exciting, but hard to review compre- 
hensively. In this paper, I have tried to touch on examples which are particularly 
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instructive and have been forced to leave out a wide range of other interesting points. 
What should be clear, however, is that the recent discoveries by neutrino experiments 
have opened up a wide range of interesting questions and opportunities. This promises 
to be a rich field of research for both theorists and experimentalists for years to come. 



I wish to thank A. Aguilar-Arevalo, G. Karagiorgi, B. Kayser, P. Nienaber, J. Spitz, 
and E. Zimmerman for their suggestions concerning this text. 
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